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We study the possibility for detecting gamma-ray emission from galaxy clusters. We consider 
1) leptophilic models of dark matter (DM) annihilation that include a Sommerfeld enhancement 
(SFE), 2) different representative benchmark models of supersymmetric DM, and 3) cosmic ray (CR) 
induced pion decay. Among all clusters/groups of a flux-limited X-ray sample, we predict Virgo, 
Fornax and M49 to be the brightest DM sources and find a particularly low CR-induced background 
for Fornax. For a minimum substructure mass given by the DM free-streaming scale, cluster halos 
maximize the substructure boost for which we find a factor of > 1000. Since regions around the virial 
radius dominate the annihilation flux of substructures, the resulting surface brightness profiles are 
almost flat. This makes it very challenging to detect this flux with imaging atmospheric Cherenkov 
telescopes since their sensitivity drops approximately linearly with radius and they typically have 
5—10 linear resolution elements across a cluster. Assuming cold dark matter with a substructure 
mass distribution down to an Earth mass and using extended Fermi upper limits, we rule out the 
leptophilic models in their present form in 28 clusters, and limit the boost from SFE in M49 and 
Fornax to be < 5. This corresponds to a limit on SFE in the Milky Way of < 3 which is too small to 
account for the increasing positron fraction with energy as seen by PAMELA and challenges the DM 
interpretation. Alternatively, if SFE is realized in Nature, this would imply a limiting substructure 
mass of Afiim > 10'' Mq — a problem for structure formation in most particle physics models. Using 
individual cluster observations, it will be challenging for Fermi to constrain our selection of DM 
benchmark models without SFE. The Fermi upper limits are, however, closing in on our predictions 
for the CR flux using an analytic model based on cosmological hydrodynamical cluster simulations. 
We limit the CR-to-thermal pressure in nearby bright galaxy clusters of the Fermi sample to < 10% 
and in Norma and Coma to < 3%. Thus we will soon start to constrain the underlying CR physics 
such as shock acceleration efficiencies or CR transport properties. 

PACS numbers: 95.35.-l-d, 95.85.Pw, 98.62.Gq, 98.65.-r, 98.70.Sa 



I. INTRODUCTION 

Dark matter (DM) has been searched for in direct de- 
tection experiments at accelerators and also in 
indirect detection experiments looking for signals in the 
cosmic-ray (CR) spectra of antiprotons, positrons, neu- 
trinos and all of the electromagnetic spectrum from radio 
waves to gamma-rays Q. So far, the improvements in di- 
rect detection sensitivity have put this method into focus, 
but the situation may change considerably the coming 
few years as the CERN LHC experiments collect data, 
and new gamma-ray detectors are being planned, such 
as the CTA Q. In fact, it has recently been pointed out 
that a dedicated ground-based gamma-ray detector 
would have potential that goes far beyond that of the 
other methods, depending on presently unknown param- 
eters in the particle physics models for DM. 

Among the astrophysical systems which will be very 



interesting to detect, and study, with gamma-ray detec- 
tors {Fermi, H.E.S.S., MAGIC VERITAS, and eventu- 
ally large detectors like CTAj3 belong galaxy clusters. 
The most promising directions in which to search for 
a gamma-ray annihilation signal (from the annihilation 
process itself, and also the accompanying bremsstrahlung 
and inverse Compton (IC) components coming from 
charged particles produced in the annihilations) are ba- 
sically three: 

1. The galactic center (g.c). This is where all nu- 
merical simulations of cold dark matter (CDM) predict 
the highest density. However, the detailed DM density 
in the very central part is difficult to predict, due a to 
a possibly very complicated interplay between baryons, 
DM, and the central galactic black hole. Also, it is a 
very crowded region with many gamma-ray sources like 
pulsars, CR-illuminated molecular clouds, and other su- 
pernova remnants, which have to be subtracted from the 
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data to extract the DM induced signal. In fact, there is 
a recent claim of an indication of a relatively light DM 

jarticle contribution to the gamma-ray flux from the g.c. 

^, but other hypotheses seem to work at least as well 

2. The dwarf spheroidal galaxies orbiting the Milky 
Way (MW); like Seguc-1, Ursa Minor, Draco, Sagittarius, 
Sculptor, Carina or Willman-1 |10Hl4l|. The problem here 
is that the nature of many of these small, dark matter- 
dominated galaxies is not entirely clear, and the velocity 
dispersion estimates are based on rather small numbers 
of stars. Confusion with star clusters and tidal disruption 
are other complications. Once a satellite dwarf galaxy is 
accreted by the MW, the outer regions are severely af- 
fected by tidal stripping. The longer a satellite has been 
part of our Galaxy, and the closer it comes to the cen- 
ter durin g it s pericentral passage, the more material is 



removed [15|. Thus the DM density profile is very un- 



certain for most of them, especially for radii larger than 
those probed by the stars. Nonetheless, by stacking the 
data together from many dwarf spheroidals these uncer- 
tainties can be made less severe, and preliminary results 
from Fermi-JjAT shows this method to give quite promis- 
ing results [l6j . 

3. Galaxy clusters. This possibility has been less stud- 
ied theoretically, however currently there is an ongoing 
observational campaign to detect gamma-ray emission 
from galaxy clusters p7l - l28t . In fact, we noted in a pre- 
vious Letter [2^ that there are certain advantages that 
work in favour of this possible target for gamma-ray de- 
tection of DM annihilation. Galaxy clusters constitute 
the most massive objects in our Universe that are form- 
ing today. This causes their DM subhalo mass function 
to be less affected by tidal stripping compared to galaxy 
sized halos that formed long ago. The annihilation lu- 
minosities of the DM halo component for e.g. the Virgo 
cluster and the Draco dwarf scales in a way (see |29j ) 
that the ratio of gamma-ray luminosities from the smooth 
components is around 4, in favour of Virgo. In addition, 
there may be a further enhancement due to substructure, 
which to a large extent should be unaffected by tidal dis- 
ruption, at least in the outer regions. According to a 
recent estimate [s^l, more massive haloes tend to have 
a larger mass fraction in subhalos. For example, cluster 
size haloes typically have 7.5 per cent of the mass within 
r2oo in substructures of fractional mass larger than 10~^, 
which is 25 per cent higher than for galactic haloesQ 

In this paper, we will investigate in detail the potential 
of several of the most promising galaxy clusters to pro- 
duce an annihilation gamma-ray yield which could be ob- 
servable with present and planned gamma-ray detectors. 



Here, we sketch out the main arguments that need to be 
considered for maximizing the expected signal-to-noise 
ratio of a promising target cluster. First, we need to max- 
imize the DM annihilation flux of an unresolved cluster, 
F = AdmJdVp'^/D^ cx M^qq/D^, where A^m depends 
on the particle physics model of DM, p is the smooth DM 
density profile, D is the luminosity distance, and a ~ 1 
if we assume universality of the DM density profile. The 
dependence of p on the halo formation epoch breaks the 
universality and slightly modifies the mass dependence, 
yielding a = 0.83 [231. Additionally, the presence of a hi- 
erarchy of substructures down to small scales [3l| and the 
potential dependence of the particle physics cross section 
on the relative DM particle velocities [l^] may further- 
more modify the scaling parameter a and shall be one of 
the focus points of this work. Second, we need to mini- 
mize the expected noise which is a sum of instrumental 
noise, galactic and cluster-intrinsic foreground. While 
the galactic foreground varies across the sky, it is typi- 
cally lower for increasing latitude (away from the galac- 
tic plane). It is thought that the gamma-ray foreground 
from clusters is dominated by CRs which are accelerated 
at cosmological formation shocks and transported over 
time although there could be a substantial contribution 
of CRs from AGNs and supernova-driven galactic winds. 
In this work we will also study (and constrain) the CR- 
induced emission from clusters to identify those objects 
that are expected to be especially dim which would im- 
ply a low CR-induccd background. Since the CR-induccd 
emission is expected to scale with the thermal X-ray 
emission of clusters (33| . we estimate gamma-ray fiuxes 
of a nearby X-ray flux complete sample of galaxy clusters 
from the Rosat all-sky survey (extended Highest X-ra 
FLUx Galaxy Cluster Sample, HIFLUGCS, [H, [H)' 
This complements and extends previous work related to 
DM in clusters and to CRs in clusters ^^E^. 

In Section 2 we discuss the theory of DM and gamma- 
rays. In particular, we focus on the leptophilic (LP) 
and supersymmetric benchmark (BM) DM models in 
this work, and outline the framework for estimating the 
gamma-ray emission from various radiative processes. In 
Section 3 we calculate the gamma-ray fiuxes and spec- 
tral distributions of four clusters that are identified as 
prime targets for DM observations (Fornax and Virgo) 
as well as to have a high CR-induced gamma-ray yield 
(Perseus and Coma). For the same clusters we derive 
the gamma-ray surface brightness profiles in Section 4. 
To extend the analysis to more clusters and increase the 
list of promising targets we estimate the fluxes from DM 
and CRs of all clusters in the extended HIFLUGCS sam- 
ple in Section 5. In Section 6 we conclude and discuss 
our results. 



^ All halo masses and length scales in this paper are scaled 
to the currently favored value of Hubble's constant, Hq = 
70km s~^Mpc~^. We define the virial mass M200 a-iid virial 
radius r2oo as the mass and radius of a sphere enclosing a mean 
density that is 200 times the critical density of the Universe pa- 



^ Note that we have added the Virgo cluster to the sample that 
we nevertheless refer to as extended HIFLUGCS catalogue in the 
following. 
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II. THEORY 

We start the section by discussing two different, but 
well motivated DM models; LP and supersymmetric DM. 
We then present the framework that is used to calculate 
the gamma-ray emission for these DM models using an 
Einasto DM density profile and the expected enhance- 
ment from DM substructures. This is followed by an 
outline of the framework of IC emission where we take 
into account photon fields of the cosmic microwave back- 
ground (CMB), dust, and starlight. We end by summa- 
rizing our formalism of calculating the CR-induced pion 
decay gamma-ray emission which is thought to dominate 
the astrophysical gamma-ray signal. 



A. Detecting Particle Dark Matter 

Besides the intrinsic interest in the gamma-ray flux 
from galaxy clusters generated by conventional hadronic 
and electromagnetic processes - which should be close 
to observability with the Fermi-LAT data [ilHiel. |47| - 
the possible contribution from WIMP DM is of greatest 
interest. A WIMP (Weakly Interactive Massive Particle) 
which fulfils the WMAP bounds on the relic density of 

CDM of ii 



^CDAih^ = 0.112 ±0.0056, 



will in many cases naturally give a gamma-yield which 
may be observable. In addition, there are possible en- 
hancement effects known, such as the astrophysical boost 
from dense substructure of DM halos, or the particle 
physics boost from the Sommerfeld effect, which may in- 
crease the chances of detection further. 

There are three methods for detecting DM candidates 
which are presently employed: First, at particle physics 
accelerators such as CERN's LHC, one is now entering 
into a new energy regime which may allow the production 
of the heavy, electrically neutral and long-lived particles 
which may constitute DM. From these experiments one 
may get the first glimpse of the mass scale beyond the 
standard model where DM may reside. However, to re- 
ally show that any of the hypothetical, new particles cre- 
ated at the LHC is actually the DM, one has to rely on the 
two other methods available for the DM search, namely 
direct and indirect detection. Second, direct detection 
methods, which are presently evolving rapidly, use the 
feeble interaction between DM particles of the galactic 
halo and nuclei such as Germanium, Sodium, Iodine, Ar- 
gon or Xenon to infer the scattering cross section and a 
rough estimate of the mass of the DM particle (for the 
currently most sensitive search, see [iol. ISOj). A charac- 
teristic of this method is that it basically only depends on 
the local DM density (which is rather well determined to 
be around 0.4 GeV/cm"^) and the basic cross section for 
DM - nucleus scattering. A disadvantage is that it does 
not benefit from the two enhancement mechanisms men- 



tioned above, boost from substructure and/or the Som- 
merfeld effect. Also, it cannot be excluded - although it 
seems at present improbable - that the local halo struc- 
ture is such that the solar system happens to be in an 
underdense region. 

Third, in indirect detection, in particular in the photon 
channel, one searches for products of DM annihilation in 
the galactic halo and beyond. Particularly interesting 
targets are the galactic center (for a recent possible in- 
dication of a signal, see [81, however, see also 0), the 
dwarf galaxies surrounding the MW [l2^51^5^. galaxy 



clusters (the focus of this work) [2^ [36 5j Sa. and even 
the cosmological large scale structure 5; 



1. Leptophilic models 

There is also a possibility to search for DM annihi- 
lation in the Milky Way indirectly through annihilation 
to antimatter channels, which however lacks the impor- 
tant directional signature of gamma-rays. Recently, it 
has however been much in focus due to the surprising 
findings of PAMELA gj] and Fermi [6^ . Viable particle 
physics models are rather constrained by other observa- 
tions, however. For example, the non-observation of an 
enhanced antiproton flux by PAMELA [6^ means that 
quark-antiquark flnal states have to be suppressed. This 
has led to the postulate of a whole class of leptophilic 
models with mainly or exclusively annihilation to leptons 
and antileptons. In addition, one has to have a sizeable 
fraction of final states containing muons or tau leptons 
as the shape of the spectrum disfavors direct annihilation 
to electrons and positrons only. A large enhancement, of 
the order of at least several hundred, of the annihila- 
tion rate is also needed, something that may be given by 
the Sommerfeld e ffect (for a pedagogical review of viable 
scenarios, see [32|)- As still another difficulty for these 
models, the lack of an IC signature from the Galactic 
center means that the DM density distribution must be 
cored rather than cuspy p37 H63l . (For a review of the DM 
modeling of these effects, see For a recent treatment, 
showing still viable models, see (70j.) 

A more standard explanation of the PAMELA/Fermz 
excess would be pulsars or other supernova remnants 
(e.g., Inl-Izl)- In that case, DM would more likely be 
explained by the conventional scenario of a WIMP, for 
example the thoroughly studie d lig htest supersymmetric 
neutralino (for reviews, see |74| - [7a |). We note that AMS- 
02 experiment on the International Space Station (77j 
may give interesting clues to the origin of the antimatter 
excess. 

In the scenario with leptophilic DM, one has to rely 
on some enhancement of the annihilation rate from the 
effects of DM halo substructure, a possibility which was 
realized long ago [TSl - lsoj . However, a boost factor as 
large as several hundred is difficult to achieve in the solar 
neighborhood, due to tidal stripping of subhalos in the in- 
ner part of the Galaxy. Sommerfeld enhancement is thus 
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important for Icptophilic models that claim to explain 
the PAMEL A/Fermi result. This effect was computed 
for electromagnetism by Arnold Sommerfeld many years 
ago [m and rece ntly rediscovered in the quantum field 
theory of DM |82l - [84j . In quantum mechanics describing 
electron scattering and electron positron annihilation, it 
is caused by the distortion of the plane wave describing 
the relative motion of the annihilating particle pair, due 
to the near formation of a bound state caused by pho- 
ton exchange. In the ladder approximation for QED, one 
reproduces the Sommerfeld effect, and the square of the 
relative wave function at the origin (which enters into the 
probability for the short-distance process of annihilation) 
is increased by the factor [131 
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(1) 



with a = e^jhc the fine-structure constant, and v the 
relative velocity. This amounts to Sqed = Tra/u for 
small velocities. For a Yukawa-like particle of mass m^, 
mediating a weak attractive force with coupling constant 
Qfy between DM particles x with mass m^, the small- 
velocity limit of the enhancement becomes instead 



Sy 



(2) 



In the general case, with mediators that may also excite 
virtual charged particles in the DM sector, one has to 
solve numerically a coupled system of differential equa- 
tions using appropriate boundary conditions {82t - l84j . In 
some cases, the enhancement factor S can be as high 
as several hundred to a few thousand, depending how- 
ever on the exact parameters of the theory. The effect is 
usually strongly velocity-dependent, depending on veloc- 
ity as 1/v or, in the fine-tuning case of being very near 
resonance, as This means that in a virialised sys- 

tem (such as a galaxy cluster) with large velocity disper- 
sion I'd the SFE will be smaller than the one expected 
in a single galaxy such as the MW, roughly by a fac- 
tor wmw/'^^cI- Note that the 1/w-scaling is valid only for 
(w/c) > (m^/m^); at smaller velocities and outside reso- 
nances, the 1/w-enhancement saturates at m^/m^ (sBj . 

As LP models, apart from being slightly contrived from 
the particle physics point of view, are also rather limited 
by several sets of astrophysical data. We take advantage 
of the recent reanalysis [70] to define a benchmark LP 
model that is still viable. It is found that the bounds from 
WMAP5 approximately imply a Sommerfeld boost factor 
S which has to satisfy S{v 150km/s) < 250/ f {m^/l 
GeV) , where / is the fraction of energy from DM annihi- 
lation that ionizes the intergalactic medium, / ~ 0.7 for 
annihilation to electrons and a factor of a few smaller for 
all other standard model final states. 

Taking into account the whole cosmic history of LP 
models, also utilizing limit on spectral and polarization 
distortions of the CMB, the maximal value for the boost 
factor for a 1 to 2 TeV particle and sub-GcV force carriers 



is found to be (T^I between 400 and 800. Although lower 
than the first estimates (e.g., [H, [s^l) this is stiU enough 
to explain the PAMELA and Fermi excess, given other 
astrophysical uncertainties. 

We adopt our benchmark LP model from [t^ where 
we use a DM mass = 1.6 TeV and a branching ratio 
of (1/4 : 1/4 : 1/2) into fi' : c+e" : tt+tt"). Further- 
more, since most of the gamma-ray flux is expected to 
come from dense subhalos within clusters which have a 
velocity dispersion close to the velocity limit Wgat where 
the SFE saturates, we use S = Bsfo(wsat) ~ 530 for the 
cluster halosQ However, for figures where we neglect the 
substructure boost, we adopt a velocity dependent SFE 
that is normalized to fit the electron and positron excess 
observed at Earth, 



Bsfo(t'cl) = Bsfo,MW— — = 70 
Vcl 



M. 



200 



1015 Mp 



-1/3 



(3) 



where the local boost factor required to explain the data 
is Bsfc,MW = 300. The velocity dispersion for the MW, 
wmw ~ 220 km s^ , and the mass dependent velocity dis- 
persion Uci for clusters is derived in [ssj . 



2. Supersymmetric dark matter 

The most studied models for DM are supersymmetric 
ones, especially models where the lightest supersymmet- 
ric particle, in most models the lightest neutralino, is 
stable. The stability is assured in viable supersymmetric 
models due to a discrete symmetry, i?-parity. This sym- 
metry is needed from the particle physics point of view 
to avoid fast proton decay, and automatically makes the 
lightest supersymmetric particle stable, and therefore a 
good candidate for DM. In addition, as the coupling to 
ordinary matter is given by gauge couplings, the neu- 
tralino is automatically a WIMP candidate. Unfortu- 
nately, the breaking of supersymmetry (which has to be 
there since otherwise, e.g., there would exist a scalar elec- 
tron with the same mass as the electron) means that a 
large number of essentially free parameters enter actual 
calculations. The neutralino is then a linear combination 
of the supersymmetric partners of the U(l) gauge bo- 
son, B, one of the SU(2) gauge bosons, W3, and the two 
Higgs doublets in the Minimal Supersymmetric extension 
of the Standard Model (the MSSM), and H^, with 
mixing coefficients which depend on the supersymmetric 
breaking parameters [89| . The parameter space is very 
large, too large in general to scan efficiently even with 
the most powerful computers. Therefore, often simpli- 
fying assumptions are made to limit the number of free 



* Note that S is strictly speaking the Sommerfeld enhanced anni- 
hilation cross section, but here we have additionally included the 
order unity constraint that ensures that the relic density is not 
over-depleted at freeze-out. 
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parameters. Even then, a large number of viable mod- 
els is found which give a relic density consistent with 
the WMAP data. This means that they are indeed very 
good templates for WIMP DM, something that explains 
their popularity besides the subjective statement that 
few more well-motivated models for DM have been put 
forward so far. 

The wave function of the lightest neutralino can be 
written 



Xi = aiB + a2W3 + 03^1° + aiH^ , 



(4) 



BM 


■mx 




(av) 






[GeV] 








I' 


140 


0.09 


4.0 X 10 


-27 


J' 


315 


0.12 


3.3 X 10 


-28 


K' 


570 


0.10 


4.4 X 10 


-26 


J* 


234 


0.09 


8.9 X 10 


-29 



TABLE I. Relevant parameters for the four benchmark mod- 
els. The mass of the DM particles are given denoted by 171^, 
the relic density by Q^h^ ^-nd the annihilation rate today by 
{av). 



with 



El 



= 1 



(5) 



where the gaugino fraction of a given neutralino is |ai | 
02!^ whereas josp -I- |a4p is the higgsino fraction. 



We use the DARKSUSYpackage [90| to compute the 
mass, couplings and relic density for a given set of pa- 
rameters. Actually, we take advantage of so-called bench- 
mark (BM) models, which have been proposed in differ- 
ent contexts (in particular in (lo| giving predictions for 
imaging air Cherenkov telescopes, like MAGIC II and 
CTA). They give a good representation of models with 
high enough gamma-ray rates to become the first candi- 
dates for DM indirect detection in imaging air Cherenkov 
telescopes. We use the following supersymmetric BM 
models: 

• I': This model was introduced in [9l|, where its 
phenomenology at colliders was studied. Annihila- 
tion directly into lepton pairs is suppressed for neu- 
tralinos due to their Majorana nature and therefore 
helicity suppression for annihilation in the galac- 
tic halo [93 • The higher order process (inner 
bremsstrahlung, IB) XiXi — ^ which does 
not suffer from helicity suppression [o^ [oj] , gives 
a considerable contribution due to light sleptons in 
this model. 

• J': This model, also from [9l| is in the so-called 
coannihilation tail. The sleptons are nearly degen- 
erate with the neutralino, which causes the large 
IB from the Icptonic final states to give a high en- 
hancement of the flux. 

• K': A representative model for the funnel region, 
where the annihilation dominantly occurs in the 
s-channel thr oug h exchange of the pseudoscalar 
Higgs boson [91|. Here IB contributions are not 
important. 

• J*: Annihilation in the so-called coannihilation re- 
gion, introduced in [o^l as BM3, with a particularly 
large IB contribution. 

The details for the DM BM models arc summarized in 
Table H 



3. Final state radiation 

There are two types of radiative processes which are 
important for our BM models, first the usual QED final 
state radiation for both the LP model and for super- 
symmetric models with charged final states. In addition, 
there may be direct emission from a virtual, charged, ex- 
changed particle (such as a spin-0 slepton or squark), in- 
ternal bremsstrahlung, IB. The latter is essentially only 
important when there is a helicity suppression for the 
lowest order annihilation process, which is the case for 
neutralinos, whose Majorana nature make the annihi- 
lation rate in the s wave proportional to for final 
state fermion /. Interestingly, both the gamma-ray and 
fermion energy spectra are peaked at the highest energy 
for this process, which can in some models cause a rather 
spectacular bump in the spectra for these models. 

The photon spectrum resulting from final state radia- 
tion is universal with only a weak dependence of the un- 
derlying particle physics model. The photon yield from 
this process is given by (see e.g. [9^ ) 



dN 



XX 



dx 



aQx TT / M 
-Tx[x) log 



Ami (1 



x) 



(6) 



Here, the normalized photon energy x ~ E^/m^c'^, 
a = /he is the fine-structure constant, Q\ and mx 
the charge and mass of the final state particle X , respec- 
tively. The function Tx{x) depends on the spin of the 
final state and is given by 



ermion 



{x) 



i + {i-xY 



(7) 



for fermions. 

The differential energy spectrum for the IB process is 
more complicated (see [93l |94| ) , but can be computed us- 
ing DarkSUSY. The differential photon energy spectra 
for all our BM models are shown in Fig. [TJ The peaking 
at high photon energy and the resulting bump caused by 
inner bremsstrahlung is clearly seen for BM models I', J', 
and J*. Despite the different particle masses and cross 
sections, the spectrum of the continuum emission shows 
a remarkably similar shape which suggests that a simple 
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rescaling of particle masses (responsible for the exponen- 
tial cutoff in the spectrum) and cross sections (spectral 
amplitude) should yield a roughly scale-invariant contin- 
uum emission spectrum except for the presence of the 
final state emission feature which depends on the de- 
tails of the specific decay channels. The right panel of 
Fig.Hshows the cumulative number of leptons (electrons 
and positrons) above a given energy for our different DM 
models which is of interested for the high-energy IC emis- 
sion of those leptons. Only LP models have energetic 
enough electrons such that the IC emission is powerful 
enough to either be constrained or detected at GeV en- 
ergies and higher. 



B. Astrophysical modeling of DM induced emission 

We now turn to the detailed modeling of the surface 
brightness profiles of DM annihilation emission and dis- 
cuss the DM profiles for the smooth distribution and sub- 
structures. 



1. General equations 

The differential photon flux within a given solid angle 
Af2 along a line-of-sight (los) is given by 



dE^ dAdtdE^ 



AO 



dV.sin^^(V',-B^), (8) 



where 
dS^ 



dE. 



(^, > E^) 



dn 



An 



dlq,^^{E^,r)A{d), (9) 



los 



and S^{ip,> E^) denotes the surface brightness above 
the photon energy E^. The integration along the line- 
of-sight /, in the direction -0 that the detector is point- 
ing, is parametrized such that the radius of the source 
r = y^7M-~D^j~~2l3^^7cos^, where Dd is the distance 
from the Earth to the center of the cluster halo and 
cos 4" = cos 6* cos V' — cos tp sin 9 smip. The angular in- 
tegration dft ~ sin dd9 dip is performed over a cone cen- 
tered around and the opening angle Afi is typically 
taken to be a few times the point spread function (PSF) 
^res- The limited angular resolution results in a prob- 
ability that a photon coming from a direction ijj^ is in- 
stead reconstructed to a direction tp, where the underly- 
ing probability distribution follow a Gaussian: 



1 



■ exp 



26*2 



where 9 = ^' — ip . 



(10) 

We denote the total differential source function by 
qsnmiE^ , f) , where we include contributions from five 
main processes; leptophilic DM annihilating through 
XX ~^ 0<A ~^ {4e or 4/z — >• 4e or 47r — > 4e} (neglect- 
ing the produced neutrinos) where the e+/c~ pairs IC 



upscatter background photons (LP-IC), leptophilic DM 
emitting final state radiation (LP-FSR), supersymmetric 
DM BM models where annihilating neutralinos generate 
e+/e~ pairs that upscatter background photons (BM- 
IC) and emit a continuum as well as final state radiation 
(BM-Cont), and CR proton induced tt" that decays into 
gamma-rays (CR-tt"). The source function is given by 



i(£'^,r) ^ qcR^^o{E^,r)+ 



iBt, 



(11) 

where the differential CR to gamma-ray source function 



is denoted by Qcr-ttO {E-^^r) (see Sec. Ill Dl and [4 If for 



further details) . The subscript i runs over the gamma-ray 
producing DM channels and the total differential boost 
factor for DM is given by: 



Bt 



i{r, Vci) 



Bsfo(wci) Bsub,i(?') for LP 
BsnhAr) forBM. 



(12) 



It is the product of enhancement factors from SFE 
Bsfo(wci) (see Eq. [3]and Sec. Ill AT|) and from substruc- 
ture enhancement over the smooth halo contribution 
Bsnb dr) ^ 1 + pLb,z(0/p'(0 El (see Eqs. MM and 
Sec. Ill B 3p . The DM source function from the smooth 
halo for each process is written in the form: 



E 

3 



dN^ 
~dE. 



-r,(r), 
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where the annihilation rate density is given by 

2 



r,(r) 



p(r) 



(13) 



(14) 



Here, the subscript j runs over all kinematically allowed 
gamma-ray producing channels, each with the spectrum 
dNjj/dEj and annihilation cross section (av)^. We de- 
note the DM mass with m 
profile with p(r). 



av)j. 

^ and the smooth DM density 



2. The smooth DM density profile 

Typically the universal Navarro-Frenk- White (NEW) 
density profile provides a good fit to both the observed 
and simulated clusters. It can be considered as a special 
case of the more general 5-paramcter profile: 



p{r) 



{r/r,f [1 + (r/rs)"]' 



S = 



7-/3 



a 



(15) 



Here, /3 denotes the inner slope, 7 is the outer slope, 
and a is the shape parameter that determines the profile 



^ Note that if the boost from substructures is S> 1, then the Som- 
merfeld enhancement approaches the constant saturated boost 
of 530. 
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FIG. 1. Source functions for different dark matter (DM) models. We show the simulated data from DarkSUSY [90| ]. 
and that generated by [9^, [9^ ; the solid lines show the fit to the data. Left panel: normalized differential continuum 
spectra for four different DM benchmark (BM) models; I' model (red circles), J' (orange squares), K' (green diamonds), 
and J* (blue triangles). We use Eq. (|BI[) to fit the continuum spectra. Right panel: number of electron and positron per 
DM annihilation above the electron energy for different DM models; I' BM model (red circles), J' BM model (orange 
squares), K' BM model (green diamonds), and J* BM model (blue triangles), leptophilic (LP) DM annihilating indirectly 
into electrons and positrons (purple +) and into muons (black x ) without considering Sommerfeld boosts. We use Eq. (|B2|) 
and Eqs. (|B3IB4|I to fit the spectra of electrons and positrons from BM and LP models, respectively. 



shape at the scaling radius = r2oo / c that characterizes 
the transition between the different power-law slopes. A 
cuspy NFW profile is given by (a,/3,7) = (1, 1,3). The 
characteristic overdensity for an NFW profile is given by 



Ps{c) 



200 pc 



log(l 



c/(l 



(16) 



where the halo mass dependent concentration parameter 
c is derived from a power-law fit to cosmological simula- 
tions with Mann > \^^^Mr., §3, 



3.56 X 



M. 



200 



1015 M, 







-0.098 



(17) 



This mass scaling agrees well with [98| for cluster-mass 
halos after converting the concentration definitions ac- 
cording to . In this work we choose to model the DM 
density by an Einasto density profile 



- 1 



a = 0.17, 



PcirSj) = P-2 exp <^ 

" V '--2 / J 

(18) 

that is slightly shallower in the center than the conven- 
tional Navarro-Frenk- White (NFW) profile, but provides 
a better fi t to recent simulated high resolution DM halos 
[loot Il0l| . It should also be noted that recent observa- 
tions of the Abell 383 galaxy cluster find a density profile 
with a shallower inner slope of /3 = 0.6 compared to an 
NFW profile, and /3 ^ 1 can be ruled out with > 95% 



confidence jl02| . These observations are based on lens- 
ing and X-ray measurements as well as the stellar velocity 
dispersion of the central galaxy. In Eq. we denote 
the density where the profile has a slope of —2 by p_2, 
and the radius by r_2 = ''aoo/c We use that the density 
p_2(c) = ps{c)/a{c) and determine a(c) through Af20o: 



dypein = M200 = 200pc 



47r ri 



(19) 



Here a = (4.16 — 4.30) for c = (3 — 10) and since a is a 
slowly increasing function with concentration c, we fix it 
for simplicity to a(c) w a « 4. 

In the recent dark matter simulation literature, it has 
become standard to characterize halos by the value, Knax, 
where the circular velocity V'circ('') = \/GM{< r)/r at- 
tains its maximum: 



K;irc(?'200) 



0.216 c 



log(l-hc)-c/(l-|-c) 



0.5 



M. 



0.32 
200 ■ 



(20) 

Especially for subhalos, this quantity seems to be more 
stable compared to the virial mass that is subject to tidal 
stripping. For comparison, we quote typical values for a 
large galaxy cluster and a galaxy group: 



A/200 = 1 X 10^5 



M200 = 4 X 10^^ Afe 



Rn 



1480 km/s, 
= 0.61r2oo, (21) 

520 km/s, 
= 0.44r2oo. (22) 
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3. Substructures 

High-resolution dissipationless DM simulations of MW 
type halos find substantial amount of substructures in 
the periphery of DM halos, while the substructures in 
the center suffer from dynamical friction and tidal effects 
depleting their central number densities. Since the DM 
annihilation rate depends on the density squared, the re- 
sulting flux from substructures is boosted compared to 
the smooth density distribution. While there is still a 
discrepancy in the literature of the exact value of the 
predicted boost factor of the DM annihilation luminos- 
ity from substructures in DM halos, this inconsistency 
starts to become resolved, appa rently con verging toward 
predictions at the high end [Sll, Il03l Il04l | . Following re- 
cent high-resolution simulations of Galaxy-sized halos, 
we adopt a b oost factor due to substructures of 230 for 
such a halo |l03| (which needs to be scaled to cluster 
halos as we will discuss below). 

The initial suggestions of a small total boos t (< i?20o) 
of order unity by the Via Lactea simulations jl05l | made 
assumptions of computing the boost assuming that the 
substructure luminosity follows that of the smooth DM 
distribution. Physically, this would imply substructure to 
follow the smooth DM distribution and a radially inde- 
pendent concentration parameter of substructures. Both 
assumptions are in conflict with tidal mass loss of satel- 
lites which is at work in simulations and cause the ra- 
dial number density to be anti-biased with respect to the 
host's mass density profile yielding a significantly flatter 
subhal o dist ribution compared to the smooth DM distri- 
bution |l04j . This increases the substructure boost pref- 
erentially in the outer parts of DM halos with a factor 
that ranges between 20 and 1000 (for dwarfs and galaxy 
clusters), hence depending sensitively on the host halo 
mass ( see F igure [21. Recent Phoenix simulations by Gao 
et al. |l06l | confirm the large enhancement due to sub- 
structures in clusters. In fact they find a boost of 1125 
for a Coma like cluster which agrees within 10% with our 
estimated boost (see e.g. Table HI)) . 

We use a double power-law function to fit the lumi- 
nosity from the smooth component of substructures (i.e. 
substructures within substructures arc not included) in- 
side radius r wh ich is determined for the Aquarius simu- 
lations [UliolEI Our best fit is given by 

ao C(Af20o) i200sm(M20o) 



Lsnb{< x) 

fix) 



(23) 



C(Af2oo) 



ai X 

0.76, ai 

M200 -A/rcs,sim 



0.95, a2 



200sim A/iim 



0.023 



-0.27, 

'M2, 



(24) 



Our approach of fitting the scahng behavior of -I/aub(< ^) di- 
rectly from numerical simulations self-consistently accounts for 
the radial depen dence of the substructure concentration due to 
tidal mass losses llOSll. 



Here, ac = 0.226, denote our fit variables, i200sm is 
the luminosity from the smooth halo without substruc- 
tures within r2oo and x — r/r20Q. Wc derive the normal- 
ization function C(Af2oo) in Eq. ([M]) from the simulations 
in reference [3l| using a value of Af200sim = 1.9 x 10^^ Mq 
for the mass of the MW halo in the simulation and 
A'jfros.sim = 1O^M0 for the mass of the smallest resolved 
subhalos in the MW simulation. The smallest mass of 
subhalos in reality, Mum, is determined from the free 
streaming length of DM at decoupling — an effect that 
erases structure on scales smaller than the free streaming 
length. In applying Eq. (|24p . we implicitly assume that 
the mass power-law scaling relation is valid down to the 
free streaming mass of DM halos. In the CDM universe, 
this is conventionally taken to be lO'^'M© [l07l . ll08l | (see 
[l09j | for a discussion of the range expected in various 
DM models). Note that potentially the power-law could 
flatten towards smaller mass scales although current sim- 
ulations show no hints of such a behavior and since Ein- 
stein's gravity is a scale-free theory, we do not expect 
such a behavior on theoretical grounds either. For DM 
halos more (less) massive than the MW we expect a larger 
(smaller) boost from substructures, simply because of the 
larger (smaller) mass range down to the minimum mass 

A/lim. 

To motivate the scaling of the substructure luminos- 
ity boost with limiting substructure mass in Eqs. (j23I 
I24p . we show how it derives from the sub structure mass 
function, diYsub/dM^ub oc Af^ub^ [3l|,[io3 and substruc- 
ture luminosity scaling with the limiting mass of satel- 
lites, Lsub oc Af['?j^. The total luminosity of substructures 
scales as 



^sub.tot 



-^sub -^sub 



-t^sub / ai\7sub-7T7 



for S = 0.674 . 



(25) 

Tidal truncation is responsible in shaping the substruc- 
ture luminosity scaling parameter S. Here we only sketch 
out qualitative arguments and leave details and com- 
parison to numerical simulations to future work. First, 
tidal effects truncate the subhalo profile primarily in the 
outer regions. As a result of this, the subhalo acquires a 
steeper logarithmic slope than the canonical 3 from the 
NFW-profile. Second, tidal stripping imposes a mass- 
dependency and should be stronger for smaller substruc- 
ture. The stripping efficiency (for an effective cross sec- 
tion a) only depends on the ambient DM density of the 
host halo, ridm, which implies a scaling of the mean free 
path of a subhalo of x,nfp oc l/(ndmc). Since there is 
less mass at any given density for smaller substructures, 
modest density inhomogeneities within the host halo seen 
by the orbiting substructures cause a stronger mass loss 
of these smaller satellites. Third, the tidal truncation 
radius should be a function of host halo radius due to 
the increasing smooth density profile for smaller radii. 
If we pick an effective host halo radius which dominates 
the luminosity contribution of the substructures and were 
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to fit the tidally truncated subhalo density profile with 
an NFW-profile, the resulting would be biased toward 
lower values implying a higher concentration compared to 
an isolated halo. In order to reproduce our value of S, we 
would need a concentration mass relation of c oc M~J^^^ 
[employing the L{M2oo, c) formula derived by reference 
[29|). As expected, this is steeper than the concentration- 
mass relation found for isolated halos on dwarf galaxy 
scales, c oc A/^°"^ |lld |. as well as galaxy and cluster 
scales, c oc MT^ oss ^ 

We now derive the squared density profile for the sub- 
structures using 



Psub(^) 



dV^ 



i200sm = ^dm / dr'47rr'^p^ (/) , 

J '"200 

1 + 02 log(a;) 



dV 



X = r/r2oo , g{x) =02+01 x" 



(26) 
(27) 



(28) 
(29) 



Here, ^dm = <lsms[E^ , R) / p{R)'^ represents the particle 



physics factor, where qsm,i{E^, R) is defined by Eq. ([131 
The different density profiles have some impact on the 
luminosity from annihilating DM, although the details 
of the density profile can be neglected compared to the 
dominating boost from substructures (assuming DM to 
be cold). In Fig. [2] we compare the radial dependence 
of the accumulative luminosity from different smooth 
cluster density profiles to the boosted luminosity due 
to substructure for different mass scales. We recalcu- 
late the overdensity, p^^ for the shallower density profile 
with /3 = 0.6 and resca le th e concentration parameter in 
Eq. ^ with 300/160 [lO^ to account for the more cen- 
tralized scale radius in cluster with a shallow inner slope. 
The emission of this profile with (3 = 0.6 is about 30% 
larger within r2oo compared to a cuspy NFW (/? = 1.0). 
This difference is built up within 0.1r2oo (i-e. close to Tg). 
Hence the slope of the central part of a cluster has little 
influence for the DM luminosity within r2oo as long as 
the degeneracy in the lensing mass measurements with 
Tg (which decreases for decreasing inner slope) has been 
taken into account. The emission from an Einasto density 
profile is about 50% larger than the cuspy NFW profile 
in the periphery of the cluster, where the difference is 
mainly built up at a few percent of ^200- The increase in 
luminosity due to substructures is negligible in the center 
of halos, but integrated out to r2oo the boost relative to 
the smooth emission profile amounts to approximately 20 
for dwarf galaxies, 200 for galaxies, and 10^ for galaxy 
clusters. We stress that these boost factors are only real- 
ized in the region around the virial radius of each respec- 
tive halo which is mostly tidally stripped for dwarfs in 
the MW. Hence a more realistic boost from substructures 
is probably much smaller in satellite dwarf galaxies. In 
addition, these boost factors are only realized for direct 
annihilation emission (continuum emission or final state 
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FIG. 2. Radial dependence of DM annihilation luminos- 
ity of smooth halo and substructures. The solid lines show 
the accumulative smooth luminosity from a cluster with the 
mass M200 = 10^^ Mq for three different density profiles; an 
Einasto profile with a = 0.17 (light blue), a cuspy NFW pro- 
file with /3 = 1.0 (thick dark blue), and a core NFW profile 
with P — 0.6 (thin purple). The dashed lines show the ac- 
cumulative luminosity from substructures for three different 
mass scales; an M200 ~ 10^"" A/q galaxy cluster (thick red), 
an M200 = 10^^ Mq galaxy (orange), and an M200 = 10* Mq 
dwarf galaxy (thin green). All luminosities have been nor- 
malized with the luminosity within r2oo from a cuspy NFW 
profile. We have assumed the standard value for the limit- 
ing substructure mass of Afum = Mq. Note the large 
expected boost from substructures in clusters (~ 1000), and 
the relative small boost in dwarf galaxies (~ 20). 



radiation) or IC scattering off homogeneous seed photon 
fields (CMB). For IC scattering of SD photons, the over- 
lap of final state leptons and SD photons is smaller which 
causes the substructure boost over the smooth emission 
to be reduced by roughly two orders of magnitudes. It 
should also be noted that the substructure boosted fiuxes 
from clusters arc much more extended than for dwarf 
galaxies, hence more difficult to detect with Cherenkov 
telescopes. In fact, it was shown in Sanchez-Conde et al. 
[ssj that inside 0.1°, where the sensitivity of Cherenkov 
telescopes is maximized, the expected DM flux including 
substructures from the brightest dwarf galaxy is about an 
order magnitude higher than the brightest cluster. How- 
ever, that work assumed a substructure boost of about 
50 in the most massive cluster halos which is more than 
an order magnitude smaller than what we use in this 
work, hence in projection we expect this difference to be 
a factor few smaller. 

In Fig. |3] we show the radial regions that domi- 
nate the DM annihilation luminosity, in particular we 
show the differential contribution to the DM luminos- 
ity per logarithmic interval in radius for three different 
mass scales. Solving for the maximum of this curve, 
d^i200sm/d(loga;)^ = in combination with Eq. ([27| . 
we find that the luminosity from the smooth NFW pro- 
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FIG. 3. Plot that shows where most of the DM annihilation 
luminosity originates. We show the differential contribution 
to the DM annihilation luminosity per logarithmic interval in 
radius which corresponds to x'^p^ /p2oo where p200 is the den- 
sity at r2oo. The solid lines represent cuspy /? = 1.0 NFW 
density profiles for three different mass scales; an M200 ~ 
10^^ Mq galaxy cluster (thick red), an M200 = 10^^ M© 
galaxy (orange), and an ^{200 ~ 10* Mq dwarf galaxy (thin 
green). The dashed lines show the contribution from sub- 
structures for the same three mass scales. All luminosities 
have been normalized with the luminosity within r2oo from 
a cuspy NFW profile. We have assumed the standard value 
for the limiting substructure mass of Miin, = 10"^ Mq. For 
the smooth profile, the majority of the flux is delivered by a 
region around rs/3 as indicated by the maximum value of the 
curves. In contrast, for substructures the emission is domi- 
nated by regions around r2oo. 



file peaks at r ~ r2oo/3c. Despite the cuspy nature of 
the density profile, the luminosity is not dominated by 
the central region but by the transition region where the 
profile steepens because of the larger volume available 
there. For large clusters with typical concentrations of 
c = r2Qo/fs — 4, the luminosity from the smooth profile 
is focused to the regime around 10% of r2oo- In contrast, 
emission from substructures is mainly contributed by the 
outer parts of DM halos. As shown in Fig. [31 the product 
of annihilation emissivity and emission volume increases 
towards r2oo and only starts to drop outside this radius. 
Note that, even though most of the substructure mass 
density has been erased in the central regions of DM ha- 
los, a cluster in projection has a significant enhancement 
due to substructures at a radius of just a few percent of 
r200- 



C. Inverse Compton emission 

In this section we outline the basics of inverse Comp- 
ton (IC) emission. As target radiation fields we consider 
cosmic microwave background (CMB) photons, and the 
light from stars and dust (SD). We derive an analytic 



model from which we can estimate the spectral and spa- 
tial distributions of SD as a simple function of cluster 
mass. 

I ■ ■ standard IC source function is given by |lllj : 



qiciE^,r) 



dVdtdE^ 



-caT 



dE, 



ph- 

2 



(£^ph) 



E, 



di?e^(£;e,r)(^G(re,,), 

(30) 

where E^ is the energy of the upscattering electrons and 
i?ph is the energy of the background photon field. We 
represent the Thomson cross section with ctt and the full 
diff eren tial Klein-Nishina (KN) cross section is captured 
by [ni: 

G(re, q)^2q\nq + il + 2q){l - q) + 1 (^09)' (1-g) ^ 

2 1 + Tcq 



where 



To = 



4£'ph-E'c 



and 



E^ 



Fc {Ec — Ej) 



(31) 



(32) 



The full KN cross section accounts for the less efficient 
energy transfer between the photon and electron once 
the energy of the Lorentz-boosted photon in the electron 



rest frame comes close to 



such that the scattering 



electron experiences a significant recoil. This results in a 
steepening of the IC gamma-ray spectrum. In the low en- 
ergy Thomson regime the IC spectrum ~ E-y 
however when Fo ^ 1 the IC spectrum steepens due to 
the KN suppression to E^"'= log{Ey). Here we denote the 
steady state electron spectrum by 

We account for two major contributions to the num- 
ber density of radiative background fields Up^; the CMB 
photons and the infra-red (IR) to ultra-violet (UV) light 
emitted by SD. The number density for the SD is given 
by Tiph = d2iVph/(dydSph)(^ph,r) = usDiEph,r)/El^ 
where the specific SD energy density usu{Eph, r) is given 
in the Appendix, Eq. (|X3)) . We model the CMB photon 
spectrum as a photon gas that is isotropically distributed 
and follows a black body spectrum with the temperature 
T = 2.73 K: 



iphl^-ph) 



d^N, 



ph 



1 



dV dEpi, 7r2 {hey exp(Sph/fcBT) - 1 



(33) 

Note that the typical energy of a black body photon be- 
fore scattering is given by (i?ph) = Cph/'^ph ~ 2.7 k^T, 
where riph and eph are the number- and energy-density 
derived by integrating nph(£'ph) and EpYinph{Eph) over 
the photon energy i^ph, respectively. 

The electrons injected from annihilating DM also suf- 
fer from diffusive and radiative losses. Hence we have to 
calculate the equilibrium spectrum of the electrons plus 
positrons denoted by drio/d-Bo in Eq. (pH)) . We derive this 
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stationary solution using the cosmic ray transport equa- 
tion (neglecting convection and re-acceleration effects) : 



d_ 

di 



d£^ 



i?(Se,x)V 



+ (?c(-Bo,x) , 



Arte 
dK 



d 
'dK 



dric 
d^ 
(34) 



where D{Ec,x.) denotes the diffusion coefficient and 
6(i?o,x) the energy loss term. The source function 
qe{Ec,x.) yields the number of electrons and positrons 
produced per unit time, energy and volume element at 
the position x: 



f 



(35) 



where the annihilation rate density T fir) is defined in 
Eq. (HH). The sum runs over the kinematically allowed 
annihilation final states /, each with a branching ratio Bf 
and a differential spectrum diV/ / dEc that represents the 
number of electrons plus positrons resulting from an an- 
nihilation event. We use the differential spectra derived 
from high-statistics simulations in [9^[9^ to compute the 
cumulative number of electrons and positrons resulting 
from neutralinos annihilate indirectly into je~ pairs 
as well as [i^ j yiT pairs in the LP model. We use Dark- 
SUSYto compute the e+/e^ spectra from our four BM 
models where only a fraction of the annihilating neutrali- 
nos is converted into electrons and positrons (see Sec. Ill Al 
and Fig. [1] for further details). 

The electrons and positrons loose their energy on a 
timescale that is shorter than the diffusive timescale in 
the ICM of galaxy clusters for cosmic ray electrons which 
is la rger than the Hubble time in our energy range (46l . 
Ill3l |. Hence, we neglect the first term of the r.h.s. in 
Eq. (|34p . and derive an expression for the equilibrium 
number density: 



dnc 



(i?c,r) = 



h{E, 



dE'^q,{E'^,r), (36) 



b{E,, r) 



' r2 

^CMB 

Stt 
4crTC 
3(771ec2)2 



g^(r) 
Stt 



El , (37) 
(38) 



Here we include the three main radiative loss processes 
for the cosmic ray electrons and positrons: (1) IC losses 
on CMB photons with the equivalent field strength of 
the CMB of ScMB = 3.24/iG (1 -f zf, where z is the 
cosmological redshift. (2) Synchrotron losses on ambient 
magnetic fields where we parametrize the magnetic field 
in the galaxy cluster by B{r) = 3/iG [no(r)/no(0)]"'^. We 
adopt a magnetic decline of ae = 0.7 in this work which 
follows from flux frozen magnetic fields. (3) IC losses on 
starlight and dust with an energy density msd('') given 
by Eq. (jA14p , where we outline the derivation in the fol- 
lowing. 

The emission of galaxy clusters at infra-red (IR) and 
ultra-violet (UV) wavelengths emerges from dust and 



starl ight in both the galaxies and the ICM (e.g. |ll4{ 
and jllSf ). Three distinctive components dominate these 
wavelengths: a central galaxy, the intra-cluster light 
(ICL), and individual cluster galaxies. We decided to 
use the accurately measured spectral shape of dust and 
starlight in the interstellar medium of our Galaxy to 
model the emission from clusters. We then normalize the 
two spectral components - far IR dust and starlight at 
wavelengths ranging from the near IR to UV - individu- 
ally by using stacked cluster data and employ a measured 
mass-to-starlight luminosity scaling relation derived from 
observations of the brightest cluster galaxy (BCG) jll6| . 

In Fig. m we characterize the spectr al sh ape through a 
flt to the galactic spectra presented in |ll4l | . The figure is 
showing the spectra at r = 0.03r2oo, which is the radius 
there the SD energy density of a galaxy cluster equals the 
energy density of the CMB black-body distribution. In- 
side this central radius the SD component is dominating, 
which is shown in Fig. [5] where we compare the energy 
densities from different radiation fields in a galaxy clus- 
ter with the mass M200 = 6.0 x 10 ^^Af©. For this figure 
we use two different profiles for the SD energy density, 
where the total profile includes the contribution from the 
ICL, the BCG and all the galaxies, while the galaxies are 
excluded in the smooth profile. To compute the IC emis- 
sion from SD, we require a non- negligible overlap of the 
relativistic Icpton distribution resulting from DM anni- 
hilation and SD. In fact, one can show with a simple 
order of magnitude calculation that the overlap, /ic-oi, 
of the photon field of individual galaxies (starlight and 
dust emission) and the smooth DM density is very small 
so that we can neglect the starlight contribution from 
galaxies to the IC emission for the reminder of this work. 



/] 



IC-ol 



A^galV^gal /light 
Klu 



< 



iVgalA/. 



gal 



10" 



(39) 



Here we assume that the exponential scale height of the 
stellar light is less than the scale radius of the galaxy halo, 
/light < {rslr2m, gai)^ ~ c'^^ ~ 10"^ and iVgaiATgai - 
O.lA/ciu- We denote the number of galaxies within r2oo 
with A''gai and the concentration of a galaxy DM halo 

with Cgal. 

We find that even for a cluster with a relative small 
mass the energy density of the SD components dominates 
over the CMB and the magnetic fields (with a central B 
field of 3 fjiG) within about 10% of r2oo- Outside this 
radius the CMB is dominating the energy density of the 
cluster. Note that while the magnetic field is always sub- 
dominant in the cluster for our assumptions, we keep its 
contribution to the total energy density for consistency. 
Also note that we extract the spatial distribution of the 
SD light in clusters from a stacked emission analysis of 
Sloan Digit al Sky Survey (SDSS) data at the redshift 
^ 0.25 llTf and do not attempt to correct for a poten- 
tial evolution in this component (c.f. Fig. W^ . 

The spatial distribution of far-IR to UV light emitted 
by SD is quite different from what is expected for the 
IC upscattered SD photons. We illustrate this in Fig. [51 
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FIG. 4. Spectral dependence of radiation fields in a cluster of 
galaxies. The black dotted line of the left peak show the spec- 
trum of CMB photons using a black body with a temperature 
of 2.73 K. The crosses of the middle and right peaks represent 
the measured spectra from stars extending from the near IR 
to UV and du st at far IR wavelengths (SD), respectively, and 
are derived in [ll4l | for a galaxy. We normalize the individual 
SD spectrum separately using the observed luminosity from 
SD in clusters. The SD luminosity is related to the cluster 
mass through Eqs. (|A5IA9[) . where we use have used the mass 
M200 = 6.0 X 1O^*M0 in this fi gure. We renormalize the SD 
spectra to the radius r = 0.03r2oo, where the smooth energy 
density of the SD light (see Fig. [5| equals the energy density 
of the CMB. The red dashed lines show the fitted SD spectral 
model. 
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FIG. 5. Spatial dependence of the energy density of radiation 
fields in a cluster of galaxies. The energy density of the CMB 
(blue dotted line) is isotropic with ticmb ~ 0.26 eVcm"'^. The 
energy density of the light from stars and dust (SD) is denoted 
by the red dashed line and the solid orange line for the total 
SD light and the smooth SD light, respectively. For compar- 
ison we show the energy density of the stars and a low dust 
model (SLD) with the solid yellow line. The SD light has been 
renormalized to a cluster with the mass M200 = 6.OxlO"M0. 
Finally we show the energy density of two magnetic field mod- 
els with a central magnetic field of 3 ^G. The magnetic 
field scales with the gas density to the power qb; dark green 
dash-dotted line (qb = 0.5) and light green dash-dotted line 
(qb = 0.7). Note that the SD radiation is dominating the 
energy density inside ~ 0.03 r2oo. 



D. Cosmic-ray induced gamma-ray emission 



where the accumulative luminosity from SD is compared 
to the IC SD photons for a M200 = 6.0 x IO^Mq clus- 
ter where the boost from substructures is excluded. The 
gamma-ray luminosity from IC upscattered SD photons 
is dominated from the inner parts where there is the 
largest overlap of the cuspy DM profile and the peaked 
SD distribution. This is in marked contrast to the accu- 
mulative SD luminosity in the optical that is dominated 
by the outer parts of the cluster. Interestingly, it rises 
as a function of radius and follows the distribution of 
an NFW mass profile outside 0.1r2oo- In the center, the 
electrons and positrons mainly cool by IC upscattering 
SD photons, hence the spatial dependence of the cooling 
cancels the distribution of the SD source function, which 
results in a distribution that approaches a density square 
profile in the center. We also find that the luminosity 
within r2oo from the total SD model is a factor three 
larger compared to the more realistic smooth SD model. 



In supernovac remnants and on scales of galaxies, there 
are convincing evidences of non-thermal populations. Es- 
pecially, in the MW, the cosmic rays are observed directly 
as well as indirectly through radio. X-ray, and gamma- 
ray emission. On larger scales of the order of few 100 kpc 
up to Mpc, there are currently a vast number of obser- 
vations of radio emission coming from radio mini halos 
in the centers of cooli ng fl ow clusters, radio relics in the 
periphery of clusters |ll8l |. as well as giant radio halos 
amounting to a total n umber of more than 50 diffuse ra- 
dio sources in clusters [ll9l . Il20l | . This type of emission 
is expected in clusters since the formation process of a 
galaxy cluster is a very energetic processes that induces 
both turbulence as well as frequently occurring merging 
and accretion shocks which both are thought to acceler- 
ate rclativistic non-thermal protons and radio emitting 
electrons to high energies. The precise origin of these 
electrons, in especially relics and giant radio halos, is 
still not settled. One possible scenario for the produc- 
tion of these electrons are hadronic CR interactions with 
ambient gas protons which results in charged and neutral 
pions that decay into electrons, neutrinos, and gamma- 



13 



10' 
10' 



t7 10° 



- 10-' 



10' 



10-^ 



; opt/IR SD total 


- y-ray IC from SD total 


- opt/IR SD smooth 


y-ray IC from SD smooth 


Mnfw, P=l 















0.01 



0.10 

^ I y onn 



1.00 



FIG. 6. Comparing the spatial dependence of emission from 
stars and dust (SD) to IC upscattered SD light. We show the 
luminosity without substructures inside radius r. The solid 
lines show the optical/IR emission from SD and the IC up- 
scattered SD light is shown by the dashed lines. The thick red 
lines show the total emission including the brightest cluster 
galaxy, the intra-cluster light and the additional intra-cluster 
galaxies, while these galaxies have been cut out in the smooth 
component shown by the thin orange lines. We normalize the 
SD light and the IC upscattered SD with the smooth luminos- 
ity within r2oo for each component. We find that the gamma- 
ray IC luminosities are dominated by the central regime, while 
the SD light is mainly build up in the outer parts of the clus- 
ter. For comparison, we show that the SD light traces the 
NFW mass profile (dotted blue) outside 0.1r2oo. 



rays (see P, m m [Til IT2ll - [l3l| Supporting evi- 
dence comes from the smoothness of the extended radio 
emission that often resembles that observed in thermal 
X-rays. This can be easily explained in the hadronic 
model since the long cooling time of cosmic ray protons 
(CRs) of order the Hubble time allows for a cluster-filling 
populati on o f CR s to build up over the formation his- 
tory m 1 1381 . Il39{ . The production of these secondaries 
depend both on the gas and CR densities in the clus- 
ter, where the CR density roughly traces the gas outside 
the core regime and is slightly enhanced in the center. 
This density scaling implies that clusters are great tar- 
gets for Cherenkov telescopes with a high sensitivity for 
the central parts of nearby clusters. Detecting the clus- 
ter gamma-ray emission is crucial in this respect as it po- 
tentially provides the unique and unambiguous evidence 
of CR populations in clusters through observing the tt" 
bump at about 100 MeV in the spectra. 



We adopt the universal spectral and spatial gamma- 
ray model developed by Pinzke & Pfrommer [4l| to es- 
timate the emission from decaying Tr^ts that dominates 
over the IC emission from primary and secondary elec- 
trons above 100 MeV in clusters. The gamma-ray for- 
malism was derived from high-resolution simulations of 
clusters of galaxies that included radiative hydrodynam- 
ics, star formation, supernova feedback, and followed CR 
physics using a novel formulation that trace the most 
important injection and loss processes self-consistently 
while accounti ng for the CR pressure in the equation 
of motion [13, Il40l Il4l| . We highlight two main un- 
certainties of the models, namely the acceleration effi- 
ciency at formation shocks and CR transport parame- 
ters. First, we note that the overall normalization of the 
CR and gamma-ray distribution scales with the maxi- 
mum acceleration efficiency at structure formation shock 
waves. Fol lowing recent observations at su pern ova rem- 
nants jl42l | as well as theoretical studies |l43l |. we as- 
sume the maximally allowed acceleration efficiency for 
strong shocks that transfers 50% of the shock-dissipated 
energy (kinetic energy corrected for adiabatic compres- 
sional heating at the shock) to CRs. These efficien- 
cies drop quickly for weaker shocks [47| which dominate 
the gravi tational energy dissipation inside galaxy clus- 
ters |l4l| . These efficiencies will have to be corrected 
down if this value is not realized at strong structure for- 
mation shocks. Second, these simulations (and by ex- 
tension the analytic model) neglect active CR transport 
such as streaming and diffusion relative to the gas, i.e., 
we assume that advective transport dominates and CRs 
are tightly coupled to the gas via magnetic fields tangled 
on sufficiently small scales which produces centrally en- 
hanced profiles. However, CR diffusion and streaming 
tends to drive the CR radial profiles towards being flat, 
with equal CR number density everywhere. While the 
CR streaming velocity is usually larger than typical ad- 
vection velocities and becomes comparable or lower than 
this only for periods with trans- and supersonic cluster 
turbulence during a cluster merger. As a consequence 
a bimodality of the CR spatial distribution may result 
with merging (relaxed) clusters showing a centrally con- 
centrated (flat) CR energy density profile jll3l |. This 
translates into a bimodality of the expected diffuse ra- 
dio and gamma-ray emission of clusters, since more cen- 
trally concentrated CR will flnd high er ta rget densities 
for hadronic CR proton interactions [ll3| . As a result 
of this, relaxed clusters could have a reduced gamma-ray 
luminosity by up to a factor of five. 



III. GAMMA-RAY SPECTRA 



An alternative scenario is the second order turbulent re- 
acceleration through the interaction of a previously injected rela- 
tivistic population of electrons by supernova driven winds or ac- 
tive galactic nuclei with plasma waves and magneto-turbulence. 



Spectrally resolved indirect DM searches have the ad- 
vantage of probing different DM models through their 
characteristic spectral distributions. To make current 
and future DM searches more effective it is important 
to know in which energy band to focus the efforts in or- 



14 



der to maximize potential DM signals over the expeeted 
background. 

We focus in this section on the spectral distribution of 
gamma-rays from clusters. In the LP model, DM anni- 
hilation radiation includes final state radiation and gives 
rise to substantial amounts of electrons and positrons 
that IC upscattcr background radiation fields to high en- 
ergies. We also consider four supersymmetric DM models 
with a high gamma-ray yield in the form of continuum 
emission and IC induced emission. In addition to the an- 
nihilating DM, we estimate the gamma-ray flux induced 
by shock-accelerated CRs. Note that all fluxes in this 
section are derived within an angle corresponding to r2oo 
and are convolved with a PSF of 0.1°. 

We compare the calculated fluxes to gamma-ray up- 
per limits (95% c.l.) set by Fermi-LAT after 18 months 
of observations [28|. In particular, to achieve a more 
reliable comparison we adopt, if nothing else is stated, 
the maximal spatially extended limits since the gamma- 
ray flux from our brightest clusters all have an angular 
extent > 1° on the sky when the boost from substruc- 
tures is included. In fact, published Fermi upper limits 
are not derived for the kind of extended emission that 
we find with our treatment of substructures, hence these 
Fermi limits may not be adequate for some of these clus- 
ters. This and other recent work highlights that the im- 
proved substructure models imply that Fermi limits will 
have to be re-calculated to accommodate for the large 
source extensions. Note, however, that for most clus- 
ters the assumption of a point source is well justified. 
The flux upper limits, that we compare to, are a func- 
tion of spectral index a. However, in the relevant energy 
range, 0.1 — 100 GeV, the spectral index varies within 
1.5 < a < 3.0. This changes the photon flux upper lim- 
its by < 50%, with Fermi-LAT being more sensitive to a 
hard spectrum j28|. 

In Fig. [7] we show the differential flux from the For- 
nax cluster which is one of the best clusters for indirect 
DM searches due to its high DM annihilation fluxes and 
low CR induced fluxes. We show the emission of four 
different supersymmetric BM models and contrast it to 
the emission induced by CRs. Comparing this emission 
to the differential gamma-ray upper limits set by Fermi- 
LAT, we find that the upper limits are not violated. The 
predicted DM flux that is dominated by the continuum 
emission from K' and I' models (shown in the left pan- 
els) are about a factor few below the upper limits, mak- 
ing it hard for Fermi to probe these kind of DM models 
in the near future without a significant improvement in 
the analysis from e.g. stacking of clusters and improved 
methods for analyzing extended sources. Furthermore, 
the gamma-ray signal induced by CRs is expected to be 
about a factor 10 below the DM continuum flux from 
the K' and I' BM models at 10 GeV. However, the IC 
emission from upscattered CMB and SD photons is at 
least a factor of a few lower for K' and a factor 1000 
lower for I', J', J*, than the expected flux from CRs above 
100 MeV, making it very hard to distinguish such a sig- 



nal from the foreground due to the similar spectral index. 
For energies below 100 MeV we expect IC emission by pri- 
mary shock accelerated electrons to be dominating |41| 
over the supersymmetric DM induced leptons. Hence in 
clusters, the IC emission from supersymmetric DM BM 
models can be neglected compared to the CR pion and 
DM continuum emission. 

For the LP DM models, however, the main contribu- 
tions to the expected extended gamma-ray flux is coming 
from the IC upscattered photons on various photon back- 
ground fields. Depending on the spatial and spectral dis- 
tribution of electrons and positrons as well as the photon 
background fields, the resulting spectral distributions of 
gamma-rays can differ greatly. Hence it is interesting to 
understand which spectral features and energy regimes 
are dominating. 

If the enhancement due to substructures in clusters 
is significant, then the distribution of electrons and 
positrons follows the radial profile of substructures out- 
side the cluster center where the smooth DM density pro- 
file is dominating. In addition, the cooling of the steady 
state electron and positron distribution is dominated by 
SD photons in the center. Especially since the SD distri- 
bution is centrally peaked and the substructure distribu- 
tion peaks in the outskirts around ^200, the overlap be- 
tween electrons/positrons and SD photons is small. This 
results in a suppression of the IC upscattered SD photons 
relative to the IC upscattered CMB photons, the final 
state radiation, and the continuum emission. However, 
if substructures are only marginally dominating over the 
smooth distribution in a cluster, the SD component be- 
comes relatively more important. 

In Fig. [8] we show the total IC emission as well as its 
individual contributions from different IC upscattered ra- 
diation fields. The left panel shows the gamma-ray emis- 
sion from the LP model and the right panel from the K' 
BM model, and for comparison we overplot the emission 
expected from the CRs. Due to the flat electron and 
positron spectrum resulting from the LP DM model and 
the smaller mean energy of the CMB compared to the 
SD, we find that the upscattered CMB photons domi- 
nate the total DM IC emission in the energy regime be- 
low 100 GeV, while the SD dominate above this energy. 
For the BM models, this transition energy is shifted to- 
wards smaller energies since the electron and positron 
spectrum has a steeper spectral distribution (see Fig.[T]). 
At the highest gamma-ray energies of about 100 GeV 
and above, the IC from starlight steepens because it 
probes the high energy tail of the electrons and more 
importantly, it suffers from the Klein-Nishina suppres- 
sion. When substructures are present, most of their flux 
resides in the outer parts of clusters. However, if we re- 
move the boost from substructures, the density profile of 
electrons and positrons is more centrally peaked, and the 
relative importance of the IC upscattered SD increases by 
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FIG. 7. Comparing the difTerential flux from different models: we show the continuum emission from DM benchmark 
(BM) models (blue dashed), electrons and positrons from DM BM models that Compton-upscatter CMB photons (green 
dotted) as well as dust and star photons (orange dash-dotted), and CR induced gamma-ray emission (red solid). Each 
panel is associated with an individual DM BM model; upper left I', lower left K', upper right J', and lower right J*. The 
emission is calculated for the Fornax cluster using a point spread function of 0.1°. The substructures boost the gamma-ray 
flux from IC upscattered CMB and continuum emission by a factor of 890 while the IC upscattered SD photons are only 
boosted by 20. Extended upper limits (< 1°) from Fermi-LAT after 18 months [l^l are also shown. In the near future 
we flnd it hard to detect even the brightest BM models, I' and K', where the continuum emission is dominating the total 
emission in the GeV energy range. 



a factor ^ 30 (see Fig. [8]lf|. A larger fraction of electrons 
(in the core) will now also cool by Compton-upscattering 
SD photons which suppresses the IC-upscattercd CMB 
hght. For comparison, we include the contribution of IC 
upscattered SD photons and to bracket the uncertainty 
in the SD model, with a smaller amount of dust. In this 
low-dust model, we have reduced the energy in dust by 
a factor 10. However, the resulting flux from the IC up- 
scattered dust in this model only decreases with a factor 



Note that for the leptophilic DM model this increase is much 
smaller since Sommcrfcld enhancement is no longer dominated 
by the low velocity DM particles living in the subhalos (sec 
Sec. Ill A II for more details). 



that is slightly smaller than 10 since the IC cooling of the 
electrons and positrons also decreases. 

Because of the large total boost factor for the LP model 
(^ 5 X 10^), we overproduce the upper gamma-ray limit 
in the 1 — 10 GeV energy interval set by fermi-LAT by 
about a factor 100. This strongly constrains the boost 
from both substructures and the SFE. Additionally, these 
constraints might improve with future more sensitive 
Cherenkov telescopes such as CTA. However, consider- 
ing that an indirect detection of DM in clusters relies on 
the boost from substructures whose main contribution 
comes from the periphery, we conclude that the wide an- 
gular extent of clusters on the sky in gamma-rays suggest 
that these sources are not ideal for Cherenkov telescopes 
since their sensitivity drops linearly with source exten- 
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It is also interesting to compare the total contribu- 
tion from the LP model, the brightest BM model (K'), 
and the CR induced emission. In Fig. IH] wc show the 
integrated flux from Fornax for our different gamma-ray 
models and compare it to the unresolved integrated flux 
upper limit on Fornax set by Fermi-hAT where they av- 
eraged the flux over the energy range 0.2 — 100 GeV as- 
suming a spectral index of 2. Again, the annihilation flux 
in the LP model is in conflict with the upper limits by 
Fermi-L AT , although only by a factor of 30 which is a fac- 
tor few less constraining than the differential flux in the 
energy range 1 — 10 GcV. As shown in this figure, the LP 
model is dominating the entire gamma-ray energy range 
up to the DM rest mass energy in this model of about 
1 TeV. Finally, we note that the flux the DM BM K' 
model is larger than the predicted emission from the CRs 
in the 0.1 — 100 GeV energy regime. Hence, for an exper- 
iment with a high sensitivity even for extended sources, 
the prospects for detecting the K' BM model over the 
expected gamma-ray background induced by CRs looks 
promising in clusters. Present day Chercnkov telescopes, 
however, have a trigger region that is smaller than the 
size of clusters, hence it has to be increased to several de- 
grees to overcome problems with background estimation. 
In addition, even though the projected CTA point source 
sensitivity (5a, 50h) shows the potential of this experi- 
ment in constraining leptophilic models as well as BM 
models with a very large neutralino mass m^^c^ > 1 TeV, 
wc flnd that analysis techniques have to be developed 
that enable the detection of extended sources without 
too much degradation of sensitivity. 

We continue by comparing the estimated differential 
flux from Fornax to three other clusters in Fig. (TUJ the 
close-by and well studied Virgo cluster, the X-ray bright 
Perseus clusters, and the massive merging Coma cluster. 
We summarize in Table |ll] the expected gamma-ray flux 
in our DM and CR models and the corresponding boost 
factors. We flnd high DM induced gamma-ray fluxes from 
the Fornax and Virgo clusters, which conflrms them as 
promising targets for indirect DM searches. Especially 
in the 1 — 10 GeV energy range it is quite striking how 
constraining the upper limits are for some of the clus- 
ters. At these energies, Fermi-LAT has its peak sensitiv- 
ity due to a combination of increasing effective area and 
decreasing source spectra as a function of energy. The up- 
per limits for Virgo and Perseus are unfortunately back- 
ground contaminated by AGN activity from M87 and 
NGC 1275, respectively, and do not gain much from the 
increased sensitivity. In addition, the angular radius of 
Virgo is about 6°, while the extended upper limit from 
Fermi-hAT is calculated assuming a cr = 1.2° radius. In 
fact, the ratio of the virial radius to the assumed ex- 
tension for the upper limits r2oo/o' for (Fornax, Coma, 
Virgo) are (3.7, 1.6, 5.3), respectively. Similarly, it is in- 
teresting to compare the radius that contains 68% of the 
flux to the assumed extension, res/a, which is given by 
(2.3, 1.0, 3.3) for (Fornax, Coma, Virgo), respectively. 



This further motivates a re-calculation of the Fermi up- 
per limits to account for the full extension of the galaxy 
clusters. We also find that the theoretical expectation of 
a high CR-induced flux in Coma and Perseus may signif- 
icantly complicate indirect DM searches in those clusters 
if the boost factors are much lower than assumed in this 
work. The Fornax cluster, however, is a great target for 
indirect DM studies because the relative low gamma-ray 
flux from CRs, absents of an active AGN, and high DM 
gamma-ray flux. 

What is the figure of merit for selecting the most 
promising cluster targets for indirect DM searches? To 
this end, we employ the luminosity-to-mass scaling rela- 
tions. The gamma-ray luminosity from the smooth den- 
sity profile is given by [29j 



CX / dVp{T) OC 



M200 c3 



[log(l + c)-c/(l + c)]2 

(40) 

and the total DM luminosity that includes boost factors 
for the LP and the BM models is given by 



i\,f 1.06 

T ~ T R rY 200 



where Bgub MSqq^ ■ 



(41) 
(42) 



We note that the IC from upscattered SD photons scales 
slightly softer with mass, hence gives rise to a larger frac- 
tion of gamma-rays in low mass clusters compared to the 
total gamma-ray flux. Also note that we have not in- 
cluded the SFE for the LP model in Eq. (gl]) since it 
saturates to a constant value in the substructures. How- 
ever, if only a small fraction of the DM resides in the sub- 
halos, i.e. the scaling of the LP DM becomes shallower 
due to the mass dependence inferred from the velocity 
dispersion: 
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(43) 
(44) 



In the lower panels of Fig. [101 we also provide the dif- 
ferential spectral index, — — dlog(dF/d£^)/dlog(i?), 
for the different emission models. We find very similar 
spectral indices for different clusters assuming a substan- 
tial boost from substructures, while for models without 
substructures the relative contribution from the upscat- 
tered dust and starlight breaks the spectral universal- 
ity for the DM models. In the 1 GcV — 1 TeV energy 
regime the CR spectral index is approximately constant 
^ (2.1 — 2.3), while the spectral index for the LP model 
varies substantially between 1.2 — 4.0, hence implying 
that gamma-ray upper limits are more sensitive to the 
specific energy regime for these models. Similarly, the 
index for the BM K' model is about 1.0 at 1 GeV and 
increases monotonically toward higher energies. To maxi- 
mize photon count statistics, experiments calculate band- 
integrated fluxes which implies a flxed spectral index over 
that energy range. Hence, we show in Table IIIII for 
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FIG. 8. Comparing the flux from different inverse Compton upscattered radiation fields. We show the differential inverse 
Compton emission induced by leptophilic DM in the left panel and by the K' benchmark model in the right panel. The 
contribution from each individual radiation field from top line to bottom line; CMB (green dashed-triple-dotted) , dust 
(purple dashed), low dust model (grey dashed), and stars (blue dashed-dotted) . The sum of the three components are 
shown with the orange solid line. The red solid lines show the CR induced gamma-ray flux. The black arrow shows the 
spatially extended differential upper limit from Fermi [2^ indicating that the LP model assumptions such as Sommerfeld 
and/or substructure boost are in conflict with the upper limit. All fluxes are calculated for the Fornax cluster within r2oo 
using a point spread function of 0.1°. For this cluster, the enhancement due to substructures from IC upscattered CMB 
and SD photons is 890. The saturated Sommerfeld boost is 530. 



TABLE II. Gamma-ray flux from various clusters within r2oo. 
Cluster K,(> lOOMeV) [phcm-^s-^]: F.,(> 100 GeV) [phcm^^ s"^]: M200''' Aum'" Bsfc'^' Beub*^' 





DM-LP(^) 


DM-BM-K'(5) 


CR-^o(«' 


DM-LP(*' 


DM-BM-K'(5' 


CR-7rO(«' 


[10" Mq] 


[Mpc] 






Coma 


8.2 X 10"* 


7.6x10"" 


4.1 X 10"° 


1.2x10"" 


2.8x 10"" 


1.5x 10"" 


12.9 


101 


530/65 


1290 


Perseus 


8.6x10"* 


7.8x 10"" 


1.5x10"* 


1.2x10"" 


2.9x10"" 


5.5x 10"" 


8.6 


79.5 


530/75 


1190 


Virgo 


1.6x 10"*^ 


1.5x10"'' 


1.5 X 10"* 


2.3x10""' 


5.4 X 10"" 


5.7x10"" 


6.9 


17.2 


530/80 


1120 


Fornax 


3.5x10-'' 


3.2x10"^" 


3.1x10"^'' 


5.1x10"" 


1.2x10"" 


1.1x10"" 


2.4 


19.8 


530/110 


890 



Notes: 

(1) T he mass of Fornax, Coma, and Perseus are derived from [s^, while the mass of Virgo is derived from 
[l44l |. The luminosity distance to Fornax, Coma, and Perseus are derived from [s^, while the distance to 
Virgo is derived from [145| |. All distances and masses assume Hq = 70 [km/s/Mpc]. 

(2) The boost due to Sommerfeld enhancement. The flrst value shows the saturated Sommerfeld boost realized 
when substructures are present, the latter is the Sommerfeld boost without substructures (see Sec. II A 1). 

(3) The boost due to substructures relative a the smooth DM distribution. We integrate the emission in a 
cylinder with a radial extent of 2.5-R200 along the line-of-sight and an angular size corresponding to r2oo (both 
measured from the cluster center). 

(4) The total gamma-ray flux from leptophilic (LP) DM where both the boost from substructures and Som- 
merfeld enhancement are included. Note that all these values are in conflict with upper limits from Fermi (see 
the following sections for detail). 

(5) The total gamma-ray flux from the supersymmetric K' benchmark (BM) model where the boost from 
substructures is included. See Sees. II A 2 and II B 3. 

(6) Gamma-ray flux induced by CR protons. See Sec. V C. 



our emission models in four different energy bands. We 
find smaller variations for the banded a.^-values, hence 
reducing the importance of the specific energy regimes. 

In this section, we have found that the LP DM mod- 
els overproduces both the spatially extended differential 
and integrated gamma-ray flux upper limits set by the 
Fermi-LAT 18 month data for several clusters. We can 



use these upper limits to constrain the boosts due to 
SFE and substructures. Furthermore, as more sensitive 
experiments and better upper limits emerge, we can start 
ruling out models that give rise to boost factors that are 
manifested in the gamma-ray flux. It is especially inter- 
esting to constrain the SFE, where a boost of the particle 
physics cross section of ~300 is required to explain the 
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ing clusters, although the current upper limits are not 
good enough to rule out the boost from substructures or 
SFE. Assuming that measured e+/e~ excess is due to 
SFE DM, then in order not to overproduce the differ- 
ential gamma-ray upper limits set by Fermi-LAT in the 
energy ranged 1 — 10 GeV, we can constrain the smallest 
mass of substructures to Miim > 10^ M©. Instead, as- 
suming the standard mass of the smallest substructures 
of Miim ~ 10"^ Mq, we can constrain the saturated boost 
from SFE to < 5 in M49 and Fornax, which corresponds 
to a maximum SFE in the MW of < 3. Hence, we con- 
clude that without substantial improvement in the mod- 
elling of these extended sources, Fermi-LAT will not be 
able to rule out LP models in their current form using 
clusters but can put impressive constraints on them. 



FIG. 9. Comparing the energy integrated flux from different 
models. We sliow the emission from CR induced emission 
(red solid), a leptophilic (LP) model that includes both final 
state radiation and IC upscattered CMB, dust and starlight 
(green dash-dotted) , and the benchmark K' model (BM) that 
includes continuum emission, and IC upscattered CMB, dust 
and starlight (blue dashed). The black arrow shows the spa- 
tially extended integrated flux upper limit set by Fermi-LAT 
again indicating challenges for the assumptions underlying the 
LP model. We also show projected CTA point source sensi- 
tivities (5cr, 50h). The emission is calculated for the Fornax 
cluster using a point spread function of 0.1°. The boost from 
Sommerfeld and substructures is about 530 and 890, respec- 
tively. 
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TABLE III. Banded gamma-ray spectral index q^.^ between 
energies Ei and £2- We show the spectral index for three 
different emission models where the boost from substructures 
is included; the CR induced emission, the leptophilic (LP) 
model that includes both final state radiation and IC upscat- 
tered CMB, dust and starlight, and the benchmark (BM) K' 
model that includes continuum emission, and IC upscattered 
CMB, dust and starlight, a^-values are derived for the Fornax 
cluster, although the variance between clusters is very small 
when we account for the enhancement from substructures. 
This suppresses the relative contribution from the upscattered 
dust and starlight component that depends on cluster mass. 



excess of electrons and positrons observed in the vicinity 
of Earth by PAMELA/fermz/H.E.S.S. In Fig.Owe use 
the estimated flux within r2oo of four bright and maxi- 
mal constraining clusters (Fornax, M49, NGC4636, and 
Coma, where we have excluded Virgo due to its large an- 
gular extent) to limit the SFE as well as the minimum 
mass of substructures where oc M2qq^^. We find 

that M49 and the Fornax cluster are the most constrain- 



IV. SURFACE BRIGHTNESS PROFILES 

The large angular extent of clusters on the sky in com- 
bination with the small PSF of most gamma-ray probes 
(^ 0.1°) suggest that we would be able to probe the dif- 
ferent spatial regimes of a cluster. Especially important 
is the spatial distribution of the gamma-ray emission as 
it biases upper limits derived for spatially extended clus- 
ters where a simple profile is often assumed. In fact, we 
will show that if DM substructures are present with the 
discussed abundances, the gamma-ray brightness profile 
is almost fiat and very different from the gamma-ray dis- 
tribution following a smooth density profile with a steep 
outer surface brightness slope of —5. 

In this section we investigate the gamma-ray surface 
brightness profiles induced by both CRs and different 
models of DM in more detail. We start by comparing 
the intrinsic surface brightness that includes substruc- 
tures from different clusters in Fig. [121 As expected, we 
find that the DM flux in all clusters follows the same uni- 
versal shape imposed by the substructures. Already at 
0.01r2oo, we are dominated by the substructures due to 
the projection of their peaked distribution the outer clus- 
ter regions (see Sec. Ill B 3l for a longer discussion). This 
implies that the spatial distribution of annihilating DM 
is independent of the type of DM model, and can not be 
used to separate different models. The surface brightness, 

Sj, of the DM scales as oc i7,smBsub-M200^^ °^ -^260°' 
which explains the factor two difference in normalization 
between the most massive cluster Coma and the least 
massive cluster Fornax with a mass ratio of about 10. 
The surface brightness induced by CRs is proportional 
to the projected squared gas profile with an additional 
enhancement in the center of cool core clusters due to 
the adiabatic compression of CRs during the formation 
of the cool core. Hence, we see a much larger variation 



We average the LP DM flux in the energy range 1 — 10 GeV 
before we compare it to the upper limits. 
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FIG. 10. Comparing the gamma-ray emission from different clusters. We show the differential flux of clusters in the upper 
panel of each figure: Fornax (upper left), Coma (upper right), Virgo (lower left), and Perseus (lower right). We show CR 
induced emission (red solid), a leptophilic (LP) model that includes both final state radiation and IC upscattered CMB, 
dust and starlight (green dash-dotted), and the benchmark K' model (BM) that includes continuum emission, and IC 
upscattered CMB, dust and starlight (blue dashed). The arrows show the spatially extended differential upper limits set 
by Fermi-LAT in the energy ranges 0.2 — 1 GeV, 1 — 10 GeV, and 10 — 100 GeV from left to right (using Gauss profiles 
for the source extension with a = (1.2, 1.0,0.8) for (Virgo, Coma, Fornax) and an X-ray-inferred King profile for Perseus) 
[28l |. We show the individual components for the BM model and LP model for Virgo and Perseus, respectively. In the 
lower panel of each figure we show differential spectral indices, aj, where dFj/dE^y ~ Ej . The fiux from the clusters is 
integrated out to r2oo using a point spread function of 0.1° and includes the boost from substructures. We find that the 
lower GeV-energy regime is most constraining due to the peak sensitivity of Fermi at these energies and use this regime 
to get upper limits on boost factors. Note that the CR-induced emission in Coma is close to the upper limits set by Fermi 
and will be tested in the upcoming years. 
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FIG. 11. Constraining boost factors using flux upper limits. 
The leptophilic (LP) gamma-ray emission is derived witiiin 
r2oo using a point spread function of 0.1°. In order not to 
overproduce the Fermi-LAT differential flux upper limits in 
the energy interval 1 — 10 GeV the boost from substructures 
and Sommerfeld enhancement (SFE) is constrained for four 
clusters; M49 (green dashed), Fornax (red solid), NGC4636 
(orange dash-dotted), and Coma (blue dotted). We indicate 
both the saturated SFE of 530 (right black arrow) as well as 
the local boost in the Milky Way of 300 (c.f. Eq. ^ that is 
required to explain the electron and positron excess observed 
at Earth with LP DM (left black arrow). If the DM interpre- 
tation is correct, we can constrain the smallest size of halos 
to be larger than 10* Mq . Contrarily, if the smallest size of 
DM halos is 10~® Mq, we can constrain the SFE to < 5 in 
M49 and Fornax. This corresponds to a maximum SFE in 
the MW of < 3, which would be too small to support the DM 
annihilation hypothesis for the PAMELA/_Fermi excess. 



in these profiles driven by the large variation of the gas 
mass fraction. It is noticeable that Fornax shows a small 
surface brightness. This is explained by the low gas den- 
sity outside the BCG and the expected small abundance 
of CRs in a small cluster in comparison to a massive clus- 
ter [m . This is the reason why Fornax is a good target 
for indirect DM searches, and superior to DM searches 
in dwarf galaxies. 

To quantify the impact of substructures on the spatial 
profiles in more detail we again turn to the Fornax clus- 
ter and show in Fig. [T3] a comparison of surface bright- 
ness profiles with and without substructures. It is quite 
remarkable how flat the profiles become when substruc- 
tures are present compared to the case without. This im- 
plies that the relative flux within the PSF in the center 
is comparable to the outer parts, which should increase 
the signal-to-noise significantly. If we assume that - for 
some reason - substructures are suppressed in the Fornax 
cluster, then the expected signal from DM is swamped by 
astrophysical backgrounds over the entire extent of the 
cluster, even though Fornax has one of the highest ratios 
of DM to CR induced gamma-ray brightness. 

In order to learn about the spatial distribution of the 



FIG. 12. Comparing the intrinsic surface brightness from 
different clusters without taking PSF effects into account. 
We show the CR induced emission (red), leptophilic emission 
(light green), and emission from the K' benchmark model 
(dark blue). The different line styles each represent a cluster; 
Fornax (solid). Coma (dotted), Virgo (dashed), and Perseus 
(dash-dotted). We include the boost from both substructures 
and Sommerfeld effect. The shape and normalization of the 
different DM models are very similar. In contrast, the CR- 
induced emission profiles have a much larger scatter due to the 
large variation of the gas fraction and expected CR fraction 
in these clusters. 
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FIG. 13. The impact of substructures on the intrinsic surface 
brightness profile o the Fornax cluster at 1 GeV. The emis- 
sion induced by CRs is denoted by the thin red solid line, 
the leptophilic model by thick light green lines, and the DM 
K' benchmark model by dark blue lines. The dashed lines 
show the brightness profiles where the boost from substruc- 
tures is included while the substructures are excluded for the 
dotted lines. The boost from Sommerfeld enhancement with 
and without substructures is 110 and 530, respectively. The 
boost due to substructures is about 890. Notice the nearly 
flat profiles when substructures are included, and the relative 
large boost at 0.01r2oo that is due to projection of peripheral 
substructures onto the cluster center. 
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DM flux in various clusters for different DM models and 
the associated gamma-ray components, we show in the 
left panel of Fig.[T4]thc brightness profiles for the smooth 
DM distribution for the same four representative clusters 
that were shown previously in Fig. [121 We find that the 
surface brightness, of the DM emission above 1 GeV 
in the outer parts of all clusters have the same shape. In 
addition, we show the spatial dependence of the individ- 
ual gamma-ray components for the Fornax cluster in the 
right panel. We see that in the BM model, the emission is 
dominated by the continuum emission. In the LP model, 
the emission at intermediate and large radii is dominated 
by IC-upscattered CMB photons while at small radii, IC- 
upscattered SD photons dominate the emission. The fact 
that the sum of both IC components resembles the pro- 
file of the BM model is not surprising: in this energy 
range, all the energy of the annihilation is imparted on 
leptons which radiate all their energy away through IC 
emission and the sum of the IC components represents a 
calorimeter for these radiating leptons. We note that the 
surface brightness due to SD upscattered photons domi- 
nates over larger central regions for energies > 1 GcV 
(cf. Fig. [5]). Furthermore, since there is no enhancement 
from substructures, the overall mass normalization of the 
LP DM model has a marginally negative trend ~ M^g^^. 
Hence, we only see a very small difference in outer parts 
of the different clusters from these type of models. The 
DM BM models scale as ^ ^^260^' hence more massive 
clusters show a slightly higher surface brightness. Finally, 
we note that the flux from DM without substructures is 
dominated by the CR induced emission for all clusters. 



In Sec. mil we saw that the spectral distribution of 
gamma-ray fiux from the LP model was dominated at 
high energies of > 100 GeV by final state radia- 
tion and IC upscattered starlight and dust, while the 
BM models were mainly dominated by the continuum 
emission. These high energies are obscrvationally im- 
portant for both Fermi-LAT and Chercnkov telescopes. 
In Fig. [151 we thus show the surface brightness above 
100 GeV for a realistic PSF of 0.1° where we include 
the boost from substructures. We investigate two clus- 
ters with an expected high DM signal; the Fornax and 
Virgo cluster. Both clusters show a very fiat brightness 
profile r < 0.1r2oo, due to both the effect from substruc- 
tures and the convolution with the PSF. The CR induced 
emission also show a smoothing of the central parts due 
to the PSF. In the outer part of the Fornax cluster the 
CR flux is suppressed compared to the DM fiux boosted 
by substructures, where the LP model gives rise to a 
gamma-ray flux that is more than three orders of mag- 
nitude larger. However, it should be noted that there 
are large uncertainties in the gas density profile which 
we use to estimate the CR induced gamma-ray surface 
brightness (see Fig. [531 for more details). 



V. POPULATION STUDIES: FLUX 
PREDICTIONS AND OBSERVATIONAL LIMITS 

In this section, we compute the expected gamma-ray 
fluxes from DM annihilation and CR interactions of the 
brightest clusters of the X-ray flux-complete sample in 
the local universe. All fluxes in this section are derived 
within an angle corresponding to r2oo; in addition we 
neglect the convolution with the PSF if nothing else is 
stated. We confront these predictions to upper limits 
obtained by Fermi 1.5-year data and conclude on the 
viability of the underlying models and perspectives for 
the next years of Fermi observations. 

A. Scaling relations 

First we focus on gamma-ray flux-cluster mass scaling 
relations for DM annihilation and CR induced emission 
in Fig. [TBI The mass scaling of the substructure boost 
steepens the intrinsically shallower DM annihilation rela- 
tion in Eq. (|4ip to oc -Afjoo^- However, the CR scaling 
relation is still considerably steeper as shown by reference 

m- 

L,{> 100 MeV) = 1.8 x 10« (^0^^) 

L,(> 10GeV) = 1.4xlO-(^^)'"phs-. 

(45) 

The CR luminosity scaling relationt0 include the IC 
gamma-ray contribution from shock-accelerated primary 
electrons as well as secondary electrons created in CR- 
proton interactions. However, the gamma-ray fiux is 
dominated by the decaying neutral pious. The differ- 
ence in the mass scaling arises from the larger contri- 
bution of primary IC emission for higher-mass clusters 
at 100 MeV that have a larger fraction of radio relics 
due to their greater mass accretion rates in comparison 
to galaxy groups — a dir ect consequence of hierarchical 
growth of structure [H, Il46l | . The steeper mass scaling 
of CR emission compared to DM annihilation already im- 
plies a general strategy to minimize the CR-induced fore- 
ground for DM annihilation and argues for very nearby 
groups. We note that these simulations do not include 
AGN feedback that is thought to furthermore redu ce th e 
baryon fraction in groups relative to that in clusters [l47| . 



The scaling relations show the flux inside r200 and do not include 
the contribution from galaxies. Also note that in this paper the 
scaling relations are normalized at 10^^ h^Q M (-,-,. In contrast, 
the corresponding scaling relations in Table 5 in |4ll | are normal- 
ized at 10-'^^ Mq instead of the mentioned 10"'^^ h:^^ Mq in 
the caption. However, all figures / tables / equations including the 
scaling relation figures in [4l| are derived for masses in units of 
10^5 /j-i Mq. 
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FIG. 14. The intrinsic surface brightness profiles without substructures at energies E > 1 GeV. We show the CR 
induced emission (red), leptophilic emission (hght green), and emission from the K' benchmark model (dark blue). Left 
panel: comparison of the surface brightness of different clusters; Fornax (solid). Coma (dotted), Virgo (dashed), and 
Perseus (dash-dotted). Right panel: comparison of different emission components in Fornax; dotted lines show the inverse 
Compton (IC) upscattered CMB photons, dashed lines show the IC upscattered photons from stars and dust, and blue 
solid lines show the continuum emission from the K' benchmark model. 




FIG. 15. Surface brightness predicted for Cherenkov telescopes at high energies. We show the emission above 100 GeV 
and include the boost from substructures. We use a point spread function of ^res = 0.1° that is typical for Cherenkov 
telescopes as well as the Fermi-LAT at this energy. Left panel shows the Fornax cluster and right panel the Virgo cluster. 
The gamma-ray emission is derived for the following components; CRs (red solid), continuum emission from the DM 
K' benchmark model (dark blue dash-dotted), as well as final state radiation (light green dashed) and inverse Compton 
upscattered dust and starlight (light green dotted) from leptophilic DM. 



The resulting smaller target density for hadronic CR in- 
teractions steepens the — AI relation of CR induced 
emission, making the case for groups even stronger. 

B. DM annihilation 

Figure [T7l compares Fermi upper limits on the gamma- 
ray flux with predictions of DM annihilation fluxes in 
the LP and BM models. We assume a boost factor due 
to substructures that has a constant contribution per 



decade in substructure mass and has a mass spectrum 
extending down to Earth masses for our LP and BM 
DM annihilation models. The gamma-ray fluxes arc cal- 
culated for those clusters where 18 months Fermi-LAT 
upper limits are derived (see [28j). These limits rule out 
the LP models in their present form with the mentioned 
assumptions in 28 clusters, and limit the boost from SFE 
to less than 5 in M49 and Fornax. Assuming universal- 
ity, this limits the Sommerfeld boost in the MW to less 
than 3 (see Fig. [TT|) . The flux level of the Fermi lim- 
its are an indirect measure of the ambient background 



23 




FIG. 16. Scaling relations of the cluster's virial mass M200 and the gamma-ray luminosity for energies E-y > 100 MeV (left) 
and > 10 GeV (right). Shown are the relations for the CR- induced emission that is dominated by pion decay resulting from 
hadronic CR interactions (solid), the leptophilic model of DM annihilations (dash-dotted) as well as the benchmark model K' 
of DM annihilation. Both DM models include the scaling of the substructure boost with cluster mass of Eq. (|4ip (dashed). 
Note that the scaling of the DM models are indicative for all DM annihilation models (as long as there is no additional mass 
scaling from, e.g., non-saturated SFE models) while the normalization depends on the particular cross section and neutralino 
mass. 



flux in the gamma-ray sky and/or the presence of strong 
point sources such as AGN in the immediate vicinity of 
the cluster position. Hence, the ratio of the Fermi upper 
limits to the predicted annihilation fluxes, -FulZ-FdM; is a 
good indication of the best cluster candidates for indirect 
DM experiments. We identify Fornax, M49, NGC4636, 
and Virgo to be prime candidates. Note however, that 
Virgo extends 12° over the sky which implies a lower sur- 
face brightness and lower signal-to-noise. Fermi limits on 
individual clusters are expected to improve as i/T/l.Syr, 
where T is the total elapsed time of the Fermi mission. 
We emphasize that the very inhomogeneous distribution 
of ^ul/^dm makes it unlikely to dramatically improve 
the limits through a likelihood stacking analysis as there 
are only a few clusters with comparably good flux ratios. 

In Figureim we show the DM annihilation flux predic- 
tions of all the clusters in the extended HIFLUGCS cat- 
alogue [15 . We find M49 and Fornax to be that bright- 
est clusters in the entire sample. While the LP model 
predicts that almost all clusters in the sample should 
be observed by Fermi after two years of data taking, 
the fluxes from the BM models are too low to be de- 
tected in the near future without improved modelling. 
Assuming the projected point source sensitivity (Scr, 
50h) of the future Cherenkov telescope array (CTA) of 
F^{> E^) = (4xl0-",2xl0-i2 2xl0-")phcm-2s-i 
at energies of = (10 GeV, 100 GeV, 1 TeV) we 
note that it will be very challenging to detect the DM 
annihilation signal from clusters without a Sommerfeld 
boost by Cherenkov telescopes. This is because the boost 
from substructures is extended while the sensitivity of 
Cherenkov telescopes scales approximately linearly with 
source extension. Hence the fluxes quoted in Fig. 1181 will 



have to be compared to a sensitivity that is scaled down 
by the ratio of cluster radius to angular resolution of 0.1° 
assuming current background subtraction techniques of 
Cherenkov telescopes. This important finding should en- 
courage the development of new methods to overcome the 
degradation of sensitivity for diffuse and very extended 
sources. Such a break-through would be needed to probe 
and potentially detect the presented class of BM models 
with a large investment of observational time. 



C. CR-induced emission 

Figu re [T^ shows flux predictions by an analytical CR 
model j4ll | for CR induced gamma-rays and compare 
those to the extended Fermz-LAT flux upper limits that 
assume King profiles which trace the X-ray emitting gas. 
For consistency, we use single- and double-beta model fits 
to ROSAT data of X-ray bright clusters for our gamma- 
ray flux estimates as an input [35| . The deprojection is 
performed following the formalism developed in e.g. jssj . 
These gamma-ray flux predictions are superior to those 
that use luminosity-mass scaling relation (a method used 
by reference psj to compare to Fermi limits) which do 
not account for the substantial scatter in the scaling re- 
lation. The analytic CR model that we employ here fol- 
lows CR-acceleration at structure formation shocks and 
their advectivc transport over cosmic time. This leads 
to a CR distribution with a universal spectral form and 
an almost universal spatial distribution. Depending on 
whether or not we account for the bias of "artificial galax- 
ies" in cosmological simulations, we derive an optimistic 
or conservative limit of the expected gamma-ray emission 
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FIG. 17. Fermi gamma-ray flux upper limits are contrasted to predicted DM gamma-ray fluxes. We show the mean 
differential flux in the energy range E-y = 1 — 10 GeV for 32 clusters. The extended Fermi-LAT upper limits are shown 
with black arrows (King proflles) and grey arrows (Gaussian proflles - Virgo, Coma, Fornax), while the point source limits 
are shown in light grey. The predicted fluxes are derived from the dominant inverse Compton-scattering of CMB photons 
in a leptophilic DM model (light red), and the continuum emission from the DM K' benchmark model (blue). Assuming 
a boost factor due to substructures that has a mass spectrum extending down to an Earth mass, the expected leptophilic 
fluxes are ruled out by upper limits in 28 of the clusters, with the strongest constraints set by M49 and Fornax. At the 
present time, we cannot constrain the benchmark models with Fermi-LAT data, although improved modelling of extended 
sources as well as stacking of clusters appears promising in testing benchmark models in the future. 



from decaying pions (bracketing our ignorance of the con- 
tribution of cluster galaxies and ram-pressure stripped 
gas parcels from infalling galaxies to the overall gamma- 
ray flux from a cluster). As discussed in Sec. Ill Dl the 
model predictions are subject to two main uncertainties, 
the acceleration efficiency at formation shocks and CR 
transport parameters and hence may be scaled down de- 
pending on specific values realized in clusters. 

The tightest Fermi limits are obtained in the 1-10 GeV 
regime due to the highest sensitivity of Fermi-hAT there. 



While the effective area of LAT increases up to these 
energies, the typical photon spectra decrease as a func- 
tion of energy: combining those effects selects this en- 
ergy range to be most sensitive. In this energy band, the 
Fermi limits close in on the conservative flux predictions 
of Virgo, Coma and Norma. Using integrated flux lim- 
its from Fermi, we arrive at very similar results. The 
insensitive limit on Perseus is due to the b right central 
source NGC1275 with an AGN in the center [l49| making 
it hard to determine a comparable limit as obtained for 
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FIG. 18. Comparing the DM annihilation flux from clusters in the extended HIFLUGCS catalogue. We show the energy 
integrated gamma-ray fluxes derived from both leptophilic DM that result in inverse Compton upscattered CMB photons 
(Hght red), and the continuum emission from the DM K' benchmark model (blue). The fluxes are calculated within 
r2oo for each of the 106 clusters included in the extended HIFLUGCS catalogue. The upper panel shows the energy 
integrated flux above 100 MeV and the lower panel that above 10 GeV, both as a function of HIFLUGCS cluster ID. 
For comparison we show the estimated point source sensitivity of the Fermi-hA[T after two years of data taking F^{> 
100 MeV) ~ 2 X IQ-^phcm-^s"^ and F^[> 10 GeV) ~ 9 x 10"" phcm'^ s'^ (dashed lines). The four brightest clusters 
are labelled yielding M49 and Fornax as the brightest targets. 



the other clusters. Additionally, recent gamma-ray ob- 
servations of the hea d-tail radi o galaxy IC310 complicate 
matters furthermore |l50l Il5l| . 

A quantity that is of great scientific interest is the CR 
pressure (Pgr) relative to the thermal pressure (Pth) as 
it may bias galaxy cluster observables, e.g., the Sunyaev- 
Zeldovich signal or hydrostatic mass estimates, and the 
resulting cosmological constraints if the analysis does not 
marginalize over a potential pressure contribution from 
CRs. Before we turn to the newly derived constraints on 
the relative CR pressure, Xcr = -PcR/^th, we study the 
Bayesian prior that different models (a simulation based 
and a simple analytic model) impose on the scaling of 
A^CR with cluster mass. 

In the simulation based analytic CR model [4l|, the 
CR pressure derives from the CR distribution function, 



f{p)dp = Cp~"dp, where C is the normalization, p = 
P/rripC is the normalized CR proton momentum, and a 
is the CR spectral index. The pressure is given by Pqr oc 
C (X C Pgas: where the volume weighted, dimensionless 
normalization scales as (C) cx M250'* for clusters with a 
mass M200 > 10" Mq [4i|. The thermal pressure Pth = 
ngask^T and is derived from the temperature in virial 

equilibrium that scales as /cbT oc GM200/ R200 oc M^qq- 
Hence, the volume weighted relative CR pressure scales 
as 



CR, sim/ 



(Pth) (fcBT) 



cx M' 



-0.22 



(46) 



which depends only weakly on cluster virial mass. This 
suggests that CRs, which end up in a cluster, are ac- 
celerated at the strongest formation shocks, i.e. during 
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FIG. 19. We contrast predictions of tlie CR induced gamma-ray emission using an analytical CR model [4j| to Fermi-IjAT 
flux upper limits. We show the mean differential flux in the energy range E-, = 1 — 10 GeV for 32 clusters. The extended 
Fermi upper limits assume King profiles and are shown with black arrows, while the point source limits are show with 
grey arrows. The blue boxes show the gamma-ray emission from CR-induced Tr^-decay, where the upper (lower) bounds 
show the estimates for an optimistic (conservative) model (see Sec. IV CI and [4l|] for details) . Note that our models are 
in perfect agreement with the derived upper limits from Fermi. Interestingly, the limits of Virgo, Norma, and Coma are 
closing in on our conservative predictions and will enforce constraints on the parameters of hadronic models such as shock 
acceleration efficiencies or CR transport properties in the coming years. 



proto-clustcr formation at high— z or in between voids 
and pre-coUapsed IGM. These CRs experience a similar 
transport history in the form of adiabatic compression 
or encounter shocks of similarly (weak) strength so that 
their distribution becomes only weakly dependent upon 
cluster mass. 

In contrast, the analytic (isobaric) CR model [s^, that 
was used by the Fermi collaboration to derive CR pres- 
sure constraints [28j . assumes Pgr = ^CR^th and con- 
strains Xcn by assuming no intrinsic cluster mass depen- 
dence of the CR pressure. While apparently very natural, 
this implies a strong cluster mass dependence of Xcr as 
the following estimate shows: 

(^CR.iso) = ^ cx \ (X Af-°-8 . (47) 



Here, we determine the mass scaling of the average gas 
density through its relation to the enclos ed g as fraction, 
/gas, and find that rigas oc /gas oc M°-^^^ [l52f . 

Figurel^shows the volume weighted relative CR pres- 
sure for 32 clusters, for which extended Fermi-hAT upper 
limits are derived (that assume King profiles). We adopt 
a constant (Xq-r) « 0.02 derived from the average of the 
large clusters in our simulation sample [4l|,|4^. Using the 
extended gamma-ray upper limits on most of the clus- 
ters in the Fermi sample, we can constrain (Xq-r) < 0.1. 
The most constraining clusters are Norma, Ophiuchus, 
and Coma with (Xcr) < 0.03. Most of the 32 clusters 
have limits on {Xcr) within a factor few from the best 
constraining clusters suggesting that stacking of clusters 
will be able to improve these limits (provided these anal- 
ysis adopt the simulation-based analytic model [4l| as a 
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FIG. 20. Limits on the relative CR pressure, {Xcn) = {PcB.)/{Pth) averaged across the cluster. Our constraints are 
obtained with the extended Fermi-LAT hmits in the energy range = 1 — 10 GeV that assume a King profile which 
matches the extension of the thermal X-ray emission (black) and for comparison with Fermi-LAT point-source limits (grey) 
[2^ . In computing those limits, we adopted our (conservative) analytic model [4^] and an averaged relative CR flux as 
obtained from our simulations. 



prior). 

In a recent work by the Fermi collaboration [28j , where 
the 1.5 year Fermi data was used to derive limits on Xcr, 
they found that the most massive clusters are on aver- 
age more constraining yielding the best constraints of 
the order of 5% (by adopting, however, point-like upper 
limits). In this work we find that medium-sized clus- 
ters provide similarly strong constraints. This difference 
is explained by the different priors on Pgr- We adopt a 
prior inferred from cosmological simulations that dynam- 
ically trace the CRs during cluster assembly suggesting 
that adiabatic transport is dominant in shaping the CR 
population in clusters and implying a weak mass depen- 
dence of XcK- Their analysis assumes a CR pressure 
that is independent of cluster mass, i.e., constraints on 
inherit the inverse mass dependence of the thermal 
pressure, Xqr oc 1/Pth oc M~°-®. While it is not clear 
which mass scaling of Xqh is realized in Nature, these 
considerations suggest that it will be critical for future 



work to account for the Bayesian prior in deriving limits 
on their pressure. 

To complete this section, we present our analytical 
model predictions for the CR-induced gamma-ray flux for 
all clusters in the extended HIFLUGCS sample in Fig. I?!] 
In the lower panel of this figure, we show the relative dif- 
ference between the gamma-ray flux computed using the 
analytical CR model to the fiux predicted by the mass- 
luminosity scaling relation. Statistically, one can view 
the analytical modeling as explicitly accounting for the 
scatter of the expected gamma-ray flux at a given mass. 
Cool core clusters have a denser core which increases the 
target density of gas protons for the hadronic reaction 
which in turn leads to a systematic increase of the ex- 
pected gamma-ray fluxes. We confirm this effect since 
all clusters with increased flux ratio of our analytical CR 
model relative to the scaling relation expectation are in 
fact cool core clusters. This suggests that the X-ray emis- 
sion should be a good proxy for the expected gamma-ray 
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FIG. 21. Comparison of the gamma-ray flux induced by CRs for the 106 clusters included in the extended HIFLUGCS cata- 
logue. The CR distribution follows the analytical model derived through hydrodynamical cosmological cluster simulations, 
where the spectral shape is independent of cluster mass [4l|] (see text for details). Gamma-ray fluxes are calculated within 
r2oo and are derived using a single or double beta profile for each cluster's gas density profile as obtained by ROSAT X-ray 
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side), both as a function of HIFLUGCS cluster ID. The eight brightest clusters are labelled. The lower panel shows the 
relative difference between the flux above 100 MeV computed using the analytical model and the gamma-ray flux predicted 
by the mass-luminosity scaling relation [ZJ]. The clusters with the largest positive offset are labelled. As expected, they 
all are cool core clusters. 



emission. Note however that there could be the counter- 
acting effect of CR streaming in relaxed clusters which 
would tend to decrease the gamma-ray luminosity [113 1. 
Future careful modelling is needed to quantify this effect. 



VI. CONCLUSIONS 

In this paper, we study the possibility for detecting 
gamma-ray emission in galaxy clusters. We consider 
benchmark (BM) and leptophilic (LP) models of super- 
symmetric dark matter (DM) as well as cosmic ray (CR) 
induced pion decay which is thought to dominate the as- 
trophysical signal from clusters. 



Dark matter annihilation. Once supcrsymmetric cold 
DM decouples from the expanding universe it streams 
and erases any potentially present smaller scales. When 
gravitational instability causes halos to collapse, this free- 
streaming scale imprints as a characteristic minimum 
halo mass Mmin into the hierarchy of structure formation 
and should still survive until the present time as a mass 
cutoff in the substructure mass function within larger 
halos. We show that the boost due to substructures 
dominates over the smooth component for halo masses 
M200 > and contributes a constant luminosity 

for each decade in subhalo mass. The ratio of halo-to- 
minimum subhalo mass M200 /-^min is maximized for the 
largest virialized systems in the Universe — galaxy clus- 
ters. The corresponding substructure boost can reach 
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values exceeding 10^ which should make nearby clusters 
the brightest DM annihilation sources after the MW cen- 
ter and hypothetical very close-by DM substructures. We 
stress that this conclusion relies on two hypotheses that 
are still open problems in particle physics and numerical 
cosmology. We assume (i) the existence of cold DM, i.e., 
the density fluctuation power spectrum extends down to 
the DM free streaming scale of about an Earth mass, 
and (ii) assume that structure formation at scales smaller 
than currently resolved ones 10^ Mq) proceeds scale- 
invariant, i.e., there is no non-linear mode coupling be- 
tween different scales. This could potentially erase struc- 
tures smaller than a characteristic one as the linear power 
spectrum approaches A^(fc) cx k^P{k) ~ const., imply- 
ing that structure on different scales may not form any 
more hierarchically, but may collapse at the same time. 

Since the mass density of substructures peaks at radii 
close to r'200, the contribution to the annihilation emis- 
sion from substructures is also dominated from these 
outer cluster regions. This implies an almost flat sur- 
face brightness profile of the annihilation emission which 
makes it necessary to have the entire cluster in the field- 
of-view to take advantage of the total substructure boost. 
This property makes it difficult to detect DM annihila- 
tion emission without an enhancement of the particle 
physics cross section over its standard value of av ^ 
3 X 10~^® cm"^ s~^ with imaging air Cherenkov telescopes 
due to the spatially extended boost from substructures 
over an angular extent of ~ 0.5 — 1° (as independently 
found by reference t38| ) . This large angular extent is not 
well matched with the sensitivity of Cherenkov telescopes 
which drops approximately linearly with radius outside 
the point spread function (with a characteristic scale of 
~ 0.1°), unless new background subtraction schemes can 
be found, e.g., for the survey mode of CTA. 

In this work, we thoroughly study the spectral emission 
characteristics of the DM emission. In general there are 
different radiation mechanism that emit in the gamma- 
ray regime, namely (i) continuum emission following the 
hadronization of the annihilating ncutralinos in BM mod- 
els which leads to the production of charged and neu- 
tral mesons which decay finally into electrons, neutri- 
nos, and gammas, (ii) inverse Compton (IC) emission by 
the leptons in the final state which up-scatter radiation 
fields from the cosmic microwave background (CMB), 
dust, and starlight, and (iii) final state radiation or inter- 
nal bremsstrahlung which dominates the highest energy 
emission close to the rest mass of the self-annihilating 
neutralinos (if present in the model under considera- 
tion)E!l 

For the first time, we compute the IC emission com- 
ponent not only from the upscattered CMB photons but 



While this final state radiation is sometimes not distinguished 
from the continuum emission contribution, we keep it separate 
as it dominates the very high-energy tail of the gamma-ray spec- 
trum in LP models while there the continuum emission is negli- 
gibly low. 



also from a realistically modeled photon field due to dust 
and starlight extending from infra-red to optical wave- 
lengths. To this end, we construct a rather simple an- 
alytic model for the spatial and spectral distribution of 
dust and starlight in clusters (which can also be used for 
the annihilation emission in galaxies after adopting a dif- 
ferent spatial profile). We find that in BM models, the 
continuum emission dominates over the IC emission of 
any available radiation field. In the LP models, where 
the continuum emission is low, the IC up-scattering of 
the CMB dominates below 60 GeV. For larger energies, 
the final state radiation dominates the emission up to 
the cutoff at m^c^ ~ 1 TeV. The IC up-scattering of the 
dust and starlight is always subdominant compared to 
the upscattered CMB (below 60 GeV) and the final state 
radiation (above 60 GeV). We stress, that the inclusion 
of the dust and starlight radiation fields are negligible 
in the cooling function for the leptons. The main reason 
why the IC emission from starlight and dust remains sub- 
dominant is the small overlap of this component with the 
peripherally peaked substructure mass density profiles. If 
for some reason, DM has less substructure, the relative 
importance of this IC component increases considerably. 

We identify Virgo, Fornax and M49 to be the best clus- 
ter/group candidates for indirect DM studies as they are 
expected to emit the brightest annihilation fluxlll More 
importantly, Fornax has a comparably low CR-induced 
gamma-ray fiux which may enable an indirect DM de- 
tection or a particularly tight limit on DM properties. 
Other nearby bright sources of the extended HIFLUGCS 
sample are Ophiuchus, and Centaurus. 

The non-detection of gamma-ray emission by Fermi 
in a total of 28 clusters/groups considerably constrains 
LP models for DM annihilation which were introduced 
to explain the increasing positron-to-lepton ratio beyond 
10 GeV. Assuming that the minimum DM substructure 
mass extends down to an Earth mass and that the DM 
annihilation flux is dominated by small scale substruc- 
ture, we limit the saturated Sommerfeld enhancement 
factor in M49 and Fornax to < 5. This corresponds to a 
value in the Milky Way (MW) of < 3, due to the larger 
Sommerfeld boost factor in smaller mass halos which 
have a lower velocity dispersion that enhances the parti- 
cle cross section (while assuming universality of the DM 
model). This would rule out LP models in their cur- 
rent form based on the non-observation of gamma-rays 
in any of the fore- mentioned clusters/groups and hence 
strongly challenge the DM interpretation of the increas- 
ing positron fraction with energy as seen by PAMELA. 
Alternatively, assuming the LP models to be correct, this 
would limit the minimum substructure mass to > 10^ Mq 
in M49 and Fornax which presents a problem f or st ruc- 
ture formation in most particle physics models |l09| . 



We excluded Virgo from our further analysis due to its large 
angular extent that needs a separate treatment. 
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Cosmic ray-induced emission. We substantially im- 
prove the modelling of the expected gamma-ray signal 
from pion-decay resulting from hadronic CR interactions 
with gas protons over previous work [28| . We employ 
the analytic model of the CR spectrum and spatial dis- 
tribution that is based on cosmological hydrodynamical 
simulations of galaxy clusters that sclf-consistently follow 
the evolution of CRs in cosmic structure formation (4]| . 
Adopting the model for all clusters in the HIFLUGCS 
sample of flux-limited X-ray clusters, we compare the ex- 
pected gamma-ray emission to that inferred from using 
an — M200 scaling relation. We note that the pre- 
dictions of the "optimistic" model in [4l| are in tension 
with Fermi upper limits for Norma, Coma, and Ophi- 
uchus. In the next 2-3 years Fermi will be able to probe 
the predictions of their "conservative" CR model. This 
will enable us to put realistic limits on a combination of 
the acceleration efficiency of CRs in structure formation 
shocks and CR transport coefhcients in clusters. 

As expected, the analytical CR model accounts for the 
"scatter" in the scaling relation and biases the gamma- 
ray flux high for prominent cool core clusters by up to a 
factor of 10 relative to the expectation from the scaling 
relation. This is due to the high target gas densities and 
associated CR densities in those cool cores that shape 
the spatial emission characteristic to be very similar to 
that observed in thermal X-rays. We caution that these 
predictions only apply for the case of negligible active 
CR transport such as CR streaming and diffusion rela- 
tive to the gas. Gamma-ray fluxes might be considerably 
smaller in those cluster that do not show an extended 
diffuse radio(-mini) halo that r nigh t point to a centrally 
concentrated CR distribution |ll3j . The CR spectrum 
E'^dN/dE shows the characteristic pion bump at ener- 



gies around 1 GeV. It is expected to dominate the DM 
annihilation signal of BM models for clusters/groups, al- 
though we find that the gamma-ray flux induced by the 
BM is higher for about 1/5 of the clusters in the HI- 
FLUGCS catalogue, where Fornax and M49 have the 
highest signal-to-noise ratio. (Note that we only include 
the cluster-intrinsic CR foreground and do not consider 
the dominant galactic and instrumental noise sources 
here.) 

Combining extended Fermi upper limits and our model 
flux predictions, we limit the relative CR pressure, 
XcR = PcR/Pth, ill 32 nearby bright galaxy clusters of 
the Fermi sample. The best limits are found in Norma 
and Coma of the order of 3%, with typical limits around 
10%. This is comparable to those inferred by the Fermi 
collaboration and mainly due to the differently as- 
sumed Bayesian prior on the CR pressure which implies a 
different cluster mass dependence of the resulting limits 
on XcR- 
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Appendix A: Detailed modelling of starlight and 
dust in clusters 



In this section, we present the details of our model of 
stars and dust (SD) in clusters of galaxies. In the first 
part, we use the measured spectral distribution of light 
from SD in a typical galaxy to derive a simple spectral 
model for SD in galaxy clusters. This is done by rcnor- 
malizing the galactic SD distribution to match the ob- 
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served luminosity from a cluster in the far infra-red (IR) 
and ultra-violet (UV) bands separately. In the second 
part, we derive the spatial distribution of the SD light 
from a stacked sample of clusters. 



1. Energy spectrum 

The cluster emission from far-IR to optical/UV is 
due to dust emission and stellar light, respectively. 
We distinguish three components: the brightest cluster 
galaxy (BCG), the intra-cluster light (ICL) and indi- 
vidual galaxies. While the latter distribution is highly 
clumped, the first two components are smoothly dis- 
tributed. We u se the far-IR to UV spectrum derived 
in reference [TT^ for a Milky Way (MW)-type galaxy to 
characterize the spectral distribution of SD in a cluster. 
Note that we only keep the spectral shape, and rcnor- 
malizc the amplitude of the spectral distribution using 
the luminosity from SD in clusters. 

In Fig. 21 we show the cosmic microwave background 
(CMB) black-body distribution together with the spec- 
tral distribution of SD light of a 6.0 x 10^^ Mq galaxy 
cluster. We fit the spectral shape of SD individually; the 
dust that peaks at about 10~^eV is fitted using a dou- 
ble power-law, while the broader spectral distribution of 
the stars that peaks at about 1 eV is fitted using a triple 
power-law. The best fit to the galactic spectrum is given 
by: 



""stars 



23 eV /1.23eV 



ph 



2.04 eV 



1.9 



ph 



0.78 eV 

Eph 



gal . 40eV ^ 0.0144 eV 



1.9 
20 



(Al) 



4.9 



0.0144 eV 

Sph 



1.9 



-I -4.9 



(A2) 



The specific energy density of the light from SD in a 
cluster can now be derived from 

usD{Eph,r) = j{r)usD{Eph) , (A3) 

where the spectral distribution of SD is given by 



UsoiEph) = i?ph 



d^iV, 



ph 



'd£;phdF 
for i = {stars, dust} . 



Y^NMoohfiEph). 



(A4) 



Here Eph denotes the energy of the SD photons, j{r) de- 
scribes the unitless spatial distribution of SD (derived in 
Sec. lA 2|) . and Ni{M2oo) denotes the mass dependent nor- 
malization for the SD. This function is determined using 



the total energy of the light from SD within r2oo, which 
is represented by E'i^vir- We relate this energy to the lu- 
minosity Li in each respective wavelength band (far-IR 
light for the dust, and optical/UV light for the stars) 



using 



C 



Jr2aa Ji ^Ph 



for i ~ {stars, dust} . 



(A5) 



Here we approximate the total energy of the photons 
within r2oo with the SD luminosity multiplied with the 
typical timescale that it takes for a photon to propagate 
through the cluster (we assume that the cluster is opti- 
cally thin). 

The luminosity from starlight is given by 



^stars(M200av) = 10 = = i ZJ~ ' ^® 



"3 



5.3 X lO^^'ergs^i 



(A6) 



where the average apparent magnitude of a cluster at z w 
0.25 in the (r-l-i)-band is given by nid ~ 15.5 |117| . This 
magnitude is derived from stacked clusters observed by 
the SDSS with the average mass M2ooav = 4.0 x lO^^A/©. 
We also use an apparent magnitude of the sun w 
—27.7, distance to the sun Dq = 4.85 x 10~^pc, lumi- 
nosity of the sun = 1.17x 10^'^ erg s~^, and distance to 
the cluster average of Dci = 1-26 x 10^ pc. Furthermore, 
to account for different cluster masses, we employ a sim- 
ple mass scaling for the luminosity from stars in Eq. (jA6p . 
Here we assume that the total starlight has the sam e hal o 
mass scaling as the brightest cluster galaxy (BCG) jll6j . 
which is a reasonable assumption since a large contri- 
bution of the smoothly distributed starlight in a cluster 
comes from the BCG. The starlight luminosity as a func- 
tion of mass is given by 



istars(.^^20o) — -^stars (A^200av ) 



200 



200av 



0.18 



(A7) 



We can now fix the unitless normalization constant for 
the stars in Eq. ()A5|) by integrating over the cluster vol- 
ume and spectral distribution of the stars: 



A'stars(-^^20o) 



/ A/; 



200 



200av 



0.18 



r2oo 6.0 X 10-9 kpc^ 
f j(r) dV 

■''•200 ■'^ ' 



(A8) 

The luminosity from dust is derived from the 
luminosity-richness scaling relation found in |l53l | 

Ldnst = 10^*-^ (%^)°'^- We pick a high richness clus- 
ter with A''2oo = 82 to normalize the luminosity since 
these clusters are less likely biased by chance coinci- 
dences. This richness corresponds to a virial mass of 
approximately A/2oodust = 6.0 x 10^^ M© jll6l |. We as- 
sume that the dust scales with halo mass in the same way 
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as the stars, and determine the normahzation eonstant 
for the dust in Eq. to 

r2oo4.0 X lO-'^kpe^ 



iVdust(M200) 



200 
200dust 



(A9) 



2. Radial profiles 



Our goal in this section is to derive a simple spatial 
model for the distribution of the light from SD in galaxy 
clusters. For this reason we use stacked cluster observa- 
tions from the Sloan Digital Sky Survey (SDSS) in the r 
and i-band that trace the starlight. For simplicity we as- 
sume that the dust traces the stars in the clusters. This 
is justified for sufficiently young stellar populations (blue 
BCG and spiral galaxies) or, in the case of the ICL, if 
the dust got stripped alongside the stars without having 
been destroyed by spallation process es th ereafter. The 
stacked surface brightness profiles in 117j | are measured 
at redshift z k, 0.25 with an average mass of the clusters 
of 4.0 X 10^^ M0. As already mentioned, three compo- 
nents contribute to the SD: diffuse ICL, galaxies, and the 
BCG in the center of clusters. However, the overlapping 
volume of the galactic light with the DM distribution of a 
cluster is very small compared to the overlapping volume 
of ICL- and BCG-hght with the DM distribution; thus, 
the relative contribution of the galaxies to the IC emis- 
sion is suppressed (c.f. Ea. l39p . In our benchmark model 
for SD, we remove the contribution from galaxies which 
corresponds to roughly 70 percent of the total SD light. 
In Fig. [22] we show the SDSS stacked brightness profiles 
as well as the fitted profilet0, where both the data and 
fits are normalized with the central surface brightness of 
the total SD component, S'gQ(O). Our spatial benchmark 
model is shown with the solid black line. 

Instead of modelling the surface brightness with a de 
Vaucouleur profile and a power-law, we use a double beta 
profile model to simplify deprojection. It is given by 

-3/3, + l/2 



'S'sd(?'±) 



2 

E 

i=l 



1 



(AlO) 

Our fit parameters for the normalized central brightness. 
Si, the core radius, rc-, and slope /3i arc given for the 
total model by 



Sl°^ = 1.0, 
5^°* = 2.3 X 10' 



1.8 kpc, 
0.19 r2oo : 



otot 
^2 



= 0.45 
0.44. 
(All) 



The mea sured brightness is converted into units of erg s ^ us- 
ing HH 5(mag"-2) = Mq + 21.6 - 2.51ogio[S{L0 pc'^)], 
where the sun 's ab solute magnitude in the r-band is given by 
Mq = 27.1 |l55f and the luminosity of the sun by Lq = 
^' ergs~^. 



10 = 



+ ICL+BCG+Galaxies 

o ICL+BCG 

. fit SD total 

. fit SD smooth 




10- 



10' 



10° 



r±/r20Q 



FIG. 22. Surface brightness profiles of the r-band {hu ~ 1 eV) 
obtained from stacked clusters in the S loan Digital Sky Survey 
(SDSS) at the redshift z ~ 0.25 pTj . The crosses show the 
total observed starlight including the diffuse intracluster light 
(ICL), galaxies, and the brightest cluster galaxy (BCG) in 
the center of the cluster. The diamonds denote the smooth 
part of the observed starlight of the ICL and the BCG. The 
solid line shows the fit to the data of the total light, while the 
smooth component is represented by the dashed line. Note 
that we use total SD brightness in the cluster center, Ss^{0), 
to normalize the SD surface brightness profiles. 



For the smooth model (our benchmark SD model), we 
find: 



SI" 



0.47, r; 



= 8.3 X 10"^ , r: 



= 3.9 kpc, 
= 0.19r2oo: 



-- 0.53 
0.54. 
(A12) 



The three dimensional spatial profile is derived by depro- 
jecting the surface brightness in Eq. (|A10[) (see e.g. [11] 
for details about the deprojection): 



2 

y 

i—1 



6ft 



1 



3ft 



S( i 3ft 



(A13) 

where B{a,b) denotes the beta-function |156| . 

The radial distribution of the energy density of SD 
governs the seed photon distribution for IC emission (to- 
gether with the CMB). The energy density from starlight 
and dust in a galaxy cluster is given by 



usu{r) 



dE, 



Ph 



ph 



dEphdV 



E, 



ph 



d£:, 



ph- 



usD(£'ph,r) 



jir) / dE, 



,gal 



200 J 



(i^ph) 



-Bph 



(A14) 

where we have used the specific energy density of the light 
from SD given by Eq. (jA3p . We compare the energy den- 
sity usj)(r) to other radiation background fields in Fig. [5] 
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for a galaxy cluster with the mass M200 = 6.0 x 10 Mq. 
We find that inside a few percent of of r2oo, the total 
energy density is dominated by the light from SD while 
the CMB is dominating elsewhere. This implies that if 
the boost from dark matter (DM) substructures is signif- 
icant, then the overlap of light from SD with the electron 
and positron distribution that trace the substructures is 
small. Hence the resulting flux from inverse Compton 
(IC) upscattered SD photons is suppressed compared to 
the IC upscattered CMB photons. 



BM model bi 62 63 64 65 65 br bg 

V 28 0.06 1.8 0.9 2.0 60 5 -1 

J' 34 0.08 5.5 0.8 2.6 250 5 -1 

K' 44 0.05 1.7 0.85 2.1 10 2 0.62 

J* 30 0.14 5.0 1.0 1.5 300 5 -1 



TABLE V. Fit parameters to the electron and positron yield 
above the electron energy from the supersymmetric bench- 
mark models. 



Appendix B: Source functions of different DM 
models 



We use DarkSUSY to derive the spectral distribution 
of the continuum emission from our four DM benchmark 
(BM) models. These spectra are shown in the left panel 
of Fig. [1] We choose to fit the source function in the 
energy regime 30MeV < E-y < m^c^ using the following 
functional form 



where E. 




4.5 



GeV 



(Bl) 



where the model specific parameters are given in Ta- 
ble lIVl 



BM model 


ffli 




a2 


0.3 


a4 




[GeV-^ cm" 










I' 


2.0 X 10 


-4 


0.36 


0.51 


701 


J' 


4.1 X 10 


-6 


0.40 


0.43 


1178 


K' 


1.7 X 10 


-4 


0.42 


0.37 




J* 


2.1 X 10 


-6 


0.34 


0.51 


490 



TABLE IV. Fit parameters to the continuum emission from 
our supersymmetric benchmark models. 

We also use DarkSUSY to derive the number of 
electrons and positrons resulting from each annihila- 
tion. These spectra are shown in the right panel of 
Fig. [TJ The electron and positron yield resulting from 
our supersymmetric BM model is fitted in the regime 



aOMeV < £^0 < TOxC using 



iVo,BM(> Se) = ^ X 



bl 
El 



1+ 



(t)' 



65 



where Eq 



E, 
GeV 



(B2) 



Fit parameters are given in Ta ble [Vl 

Finally, we use the data from [95|,[9g to derive the elec- 
tron and positron yield resulting from leptophilic (LP) 
DM annihilating indirectly into decaying muons or elec- 
trons and positrons. In the energy regime 30 MeV < 
Ec < m^<? we find the following best fit for the annihi- 
lation into muons: 

A^c,lpm(> -Ec) = exp [ci +C2X + cz[2- xf^] , 
E, 



X 



where 
ci = 1.49, C2 



and 



-5.00, C3 = -2.46, C4 



-4.40 . 



(B3) 



Similarly, for the annihilation into electrons and 
positrons we find the following best fit: 

,,2\ 



^c,LPc(> -Be) ^ 

where y 



exp (ci + C2 y -f C3 y 
log (1 - x) , 



Eo, 



and 



ci = 1.40, C2 = 1.11, C3 = 0.058. 



(B4) 



The LP model, that we are using in this work has a 
branching factor branching ratio of (1/4 : 1/4 : 1/2) into 
(/i+/i~ : e+e~ : tt+tt"). The contribution from the elec- 
trons from the decaying pious is relative small, hence the 
resulting electron and positron yield for the LP model is 
approximately given by 



iVcXp(> E,) = A^cXPc(> E,)/^ + iVc,LPM(> ^c)/4 . 



(B5) 



Appendix C: Gas density profiles 



The production of gamma-rays and secondaries from 
hadronic CRs in clusters depend both on the gas density 
and the CR number density, which roughly trace the gas 
outside the core regime and is slightly enhanced in the 
center. These gas density profile of galaxy clusters are de- 
rived from X-ray observations, where the measurements 
are mainly sensitive to the central parts of the clusters 
where there X-ray flux is high, while the core and outer 
parts are more difficult to measure due to the low signal 
to noise and often involves extrapolation from the cen- 
ter. The uncertainties in measurements and modelling, 
in combination with different specifications of active X- 
ray satellites, give rise to important differences in the gas 
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\w^X LAO vC/i. 




Pgas = |E,nf(0) 




N 1/2 

-3/3,- ^ -^'^ 






R pfprpri pp 




ni(0) [el/cm^; 


rc,i [kpc] /3i 722(0) [el /cm^] rc,2 [kpc] 


^2 n3(0)[el/cm3 


rc,3 [kpc] 






Fornax 


0.35 


0.36 0.54 2.2 x 10"^ 


190 


41 5.4 X 10"" 


183 


0.8 


[157] 


Coma 


3.55 X 10"^ 


245 0.65 










[35j 


Virgo 


0.15 


1.6 0.42 1.2 X 10"^ 


20 


0.47 






[158] 


Perseus 


46 X 10" 


'='[l + (r/57kpc)2]-i-«+4.8 


X 10"^ [l + (r/200kpc)2]~'' *^ 


el /cm=*] 




[159] 



TABLE VI. Electron number density profiles for selected clusters. 



profiles. Since the gamma-ray flux from a galaxy cluster 
is sensitive to the details of the gas profile, we choose to 
model the density profile of four bright clusters (Fornax, 
Virgo, Coma, Perseus) in more detail where we use the 
most recent and detailed modelled gas profile available 
in the literature. We show these profiles in table IVll 

The X-ray emitting gas in Fornax does not follow a 
simple /3-profile. Instead, based on deep Rosat data and 
supported by Chandra data, it is best modeled by a mul- 
ticomponent bidimcnsional model |l57l |. consisting of: (1) 
a central component (r < 5 kpc); (2) a "galactic" com- 
ponent (5 kpc < r < 40 kpc); and (3) an elliptical ICM 
component (r > 40 kpc) . 

In Fig. [23] we show the data points for the electron 
number density (rio) in Fornax, d erive d from the deep 
Rosat data presented in reference 157{ . and the best fit 



density profile together with the individual components. 
The total ric profile is derived from their fitted central and 
galactic surface brightness components while we re-fit the 
ICM component. The reason for the re-fitted outer com- 
ponent is the large uncertainty in the data points outside 
(0.2 — 0.3)^200 which we exclude in our fit. Instead we 
assume that the outer slop of tIq follows the outer slope 
of Fornax in the HIFLUGCS catalogue [s^. Deproject- 
ing the fitted surface brightness components yields the 
following electron number density profile: 



nc{r) 



-3ft 



1/2 



where 

ni = 0.35 clcm^'^ , 

n2 = 2.2 X lO^^elcm"^ , 

713 = 5.4 X lO^^elcm"^ , 



rc,i = 0.36 kpc, Pi 
r,,2 = 190kpc, /32 
rc,3 = 183kpc, /33 



(CI) 



: 0.54, 

41, 
0.8. 

(C2) 



Note that we remove the contribution of the flat central 
component of the fit (i = 1) in the outer parts of the 
cluster, thus we neglect the central contribution for r > 
40 kpc. 

For 106 X-ray bright clusters extended HIFLUGCS 
catalogue we adopt more general density profiles derived 
in the recent paper by Chen et al. |35| . They provide 
single beta-model density fits for all clusters, as well as 



_ 10° 



c 
o 




Fornax: 

— Hybrid profile; total 

Hybrid profile; components 

- - Rosat/HIFLUGCS profile 
Deep Rosat; Paulillo et al. 2001 



r/ r„ 



FIG. 23. Comparing different electron number density profiles 
for the Fornax cluster. Black crosses show the density profile 
inferred from dep rojec ted deep Rosat X-ray surface bright- 
ness observations [l57| . The total hybrid profile shown by 
the red solid line represent the best fit to the data, where the 
fitted individual density components are shown by the black 
dotted lines. The blue dashed line shows the single beta den- 
sity profile inferred from the HIFLUGCS catalogue. Due to 
insufficient sensitivity of the data to the outer cluster part, 
we use the outer slope from the HIFLUGCS catalogue [351 ]. 



the double beta-model fits to the surface profiles. We fol- 
low the dcprojection procedure in [ssj . and only choose 
the fitted double beta profiles with a that is smaller 
compared to the single beta profile. [^ 



Appendix D: Flux tables for the HIFLUGCS 
catalogue 

In this section we present the gamma-ray flux predic- 
tions using the clusters in the extended HIFLUGCS cat- 
alogue. We show the brightest 50 clusters in descending 



We adopt the single beta density profile for A2255 since it har- 
bors a giant radio halo and most likely has a disturbed mor- 
phology that might have given the fit the impression of a cool 
core. Furthermore, the difference in between the single- and 
double-beta profiles is marginal for this cluster. 
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order; CR proton induced tt that decays into gamma- 
rays m Table rvm the supersymmetric DM BM-K' model 
where the Neutralino emit a continuum as well as final 
state radiation in Table IVlIIl and LP DM that emit final 
state radiation and annihilate either indirectly to /e~ 
or through [i^ j [i" to c+/c^ that IC upscatter CMB pho- 
tons in Table IIXI The flux from DM has been boosted 



by the substructures, where we in addition include the 
Sommerfeld enhancement for the LP model. For com- 
pleteness, we show in Table |X] the parameters that we 
use to derive the fluxes for the clusters in the HIFLUGCS 
sample. 
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TABLE VII. Gamma-ray CR-vr' flux within r2oo from brightest 50 clusters in HIFLUGCS catalogue. 



Cluster 




F (2) 


F (3) 


F (4) 


E^dF^/dE./'''' 










(> 100 MeV) (> IGeV) (> 100 GeV) 


(5 GeV) 


(5 GeV) 


(5 GeV) 


PERSEUS 


68 


14.28 


19.15 


52.49 


16.26 


16.25 


10.06 


OPHIUCHU 


96 


7.71 


10.34 


28.33 


8.78 


8.78 


6.01 


COMA 


34 


4.15 


5.56 


15.24 


4.72 


4.72 


2.55 


NORMA 


94 


2.98 


3.99 


10.94 


3.39 


3.38 


2.19 


CENTAURU 


30 


2.61 


3.50 


9.58 


2.97 


2.97 


2.62 


TRIANGUL 


95 


1.77 


2.37 


6.50 


2.01 


2.01 


0.37 


AWM7 


67 


1.75 


2.35 


6.43 


1.99 


1.99 


1.24 


A2319 


98 


1.71 


2.30 


6.30 


1.95 


1.94 


0.28 


A3571 


39 


1.61 


2.15 


5.90 


1.83 


1.82 


0.41 


3C129 


70 


1.35 


1.82 


4.98 


1.54 


1.54 


0.71 


A2199 


51 


1.11 


1.49 


4.10 


1.27 


1.27 


0.33 


2A0335 


11 


1.04 


1.40 


3.83 


1.19 


1.19 


0.19 


A2029 


43 


1.04 


1.39 


3.81 


1.18 


1.17 


0.08 


A0496 


18 


1.02 


1.36 


3.73 


1.16 


1.15 


0.27 


HYDRA 


25 


0.95 


1.27 


3.48 


1.08 


1.08 


0.84 


A0085 


1 


0.93 


1.25 


3.44 


1.06 


1.06 


0.12 


A2142 


48 


0.92 


1.24 


3.39 


1.05 


1.04 


0.06 


A3266 


17 


0.92 


1.23 


3.37 


1.05 


1.04 


0.16 


A1795 


40 


0.91 


1.22 


3.35 


1.04 


1.04 


0.11 


PKS0745 


77 


0.82 


1.10 


3.02 


0.94 


0.88 


0.03 


A2256 


54 


0.78 


1.04 


2.86 


0.89 


0.88 


0.11 


A3558 


37 


0.75 


1.01 


2.77 


0.86 


0.86 


0.12 


A3667 


56 


0.71 


0.96 


2.62 


0.81 


0.80 


0.07 


A1367 


26 


0.70 


0.94 


2.59 


0.80 


0.80 


0.43 


A0478 


14 


0.68 


0.91 


2.48 


0.77 


0.75 


0.03 


A0401 


8 


0.61 


0.82 


2.25 


0.70 


0.68 


0.05 


A0262 


5 


0.57 


0.77 


2.10 


0.65 


0.65 


0.26 


HYDRA-A 


24 


0.56 


0.75 


2.05 


0.63 


0.63 


0.05 


A0754 


23 


0.52 


0.70 


1.92 


0.59 


0.59 


0.10 


A4038 


62 


0.52 


0.69 


1.90 


0.59 


0.59 


0.13 


A2052 


44 


0.51 


0.68 


1.87 


0.58 


0.58 


0.09 


A2147 


49 


0.50 


0.66 


1.82 


0.56 


0.55 


0.11 


A0119 


2 


0.48 


0.64 


1.76 


0.54 


0.54 


0.10 


A0644 


78 


0.47 


0.63 


1.73 


0.54 


0.53 


0.04 


A1644 


31 


0.43 


0.58 


1.60 


0.49 


0.49 


0.07 


A3158 


13 


0.39 


0.53 


1.45 


0.45 


0.44 


0.04 


FORNAX 


10 


0.36 


0.49 


1.34 


0.42 


0.42 


0.41 


A2063 


47 


0.36 


0.48 


1.33 


0.41 


0.41 


0.06 


MKW3S 


45 


0.35 


0.47 


1.29 


0.40 


0.40 


0.04 


A3112 


9 


0.34 


0.46 


1.26 


0.39 


0.38 


0.02 


ZwC11742 


97 


0.34 


0.46 


1.25 


0.39 


0.38 


0.03 


A4059 


63 


0.33 


0.45 


1.23 


0.38 


0.38 


0.05 


A1651 


33 


0.31 


0.41 


1.13 


0.35 


0.34 


0.02 


A0399 


7 


0.31 


0.41 


1.13 


0.35 


0.34 


0.02 


NGC1550 


15 


0.30 


0.40 


1.09 


0.34 


0.34 


0.16 


ANTLIA 


79 


0.29 


0.39 


1.06 


0.33 


0.33 


0.22 


A2657 


61 


0.29 


0.38 


1.05 


0.32 


0.32 


0.06 


A0539 


71 


0.28 


0.38 


1.03 


0.32 


0.32 


0.07 


A3581 


41 


0.28 


0.37 


1.03 


0.32 


0.32 


0.07 


A2204 


52 


0.28 


0.37 


1.02 


0.32 


0.21 


0.01 



Notes: (1) The HIFLUGCS cluster ID. (2) In units of lO^^ph cm"^ s"\ (3) In units of 10""'ph cm"^ s"\ (4) 
In units of 10~^'^phcm~'^ s~^. (5) In units of lO^^'^ergcm"^ s~^. (6) The flux within 0.1°, smoothed with a 
point spread function of 0.1°. (7) The flux within 1.0°, smoothed with a point spread function of 0.1°. 
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TABLE VIII. Gamma-ray BM-K' Continuum flux within r2oo from brightest 50 clusters in HIFLUGCS catalogue. 



Cluster 


ID(i) 


F (2) 


F (2) 


p (3) 




£^dF.,,o.i/dS^(*'5) 


di?,,,i.o/dS^(^'6) 


S/N(' 






(> 100 MeV) (> IGeV) (> 100 GeV) 


(5 GeV) 


(5 GeV) 


(5 GeV) 




FORNAX 


10 


31.88 


on oi 
zU.zl 


11 1 

11. lo 


32.00 


31.98 


30.27 


9.93 


M49 


81 


1 n /I Q 
19.48 


1 o o 
Iz.oO 


o 1 
0.81 


19.55 


19.54 


18.13 


9.29 


OPHIUCHU 


96 


1 C 1 Q 

18. io 


11 A<\ 

11. 4y 


f; o /I 
O.o4 


18.19 


18.14 


12.31 


1.66 


CENTAURU 


30 


1 A on 

14. yy 


n [^n 

y.ou 


c; 0/1 
0.z4 


15.05 


15.02 


12.40 


2.32 


HYDRA 


25 


9.19 


K OO 

o.ez 


o oi 

o.zl 


9.22 


9.20 


6.89 


2.32 


NGC4636 


29 


y.iz 


c; 7Q 

0. to 


o 1 n 
0.19 


9.15 


9.14 


8.19 


5.38 


NORMA 


94 


o.lo 


1 7 
0.1 ( 


o o 
z.eO 


8.18 


8.16 


5.42 


1.20 


PERSEUS 


68 


( .oU 


4.90 


O 70 
Z. < O 


7.83 


7.80 


4.95 


0.53 


COMA 


34 


( .08 


A Q 1 

4.81 


Z.OO 


7.61 


7.58 


4.41 


0.95 


AWM7 


67 


( .ol 


4.04 


z.OD 


7.34 


7.31 


4.66 


1.39 


A1367 


26 


D.D / 


A oo 

4.ZO 


o oo 
Z.OO 


6.69 


6.66 


3.84 


1.95 


NGC5813 


91 


D./O 




o 1 n 
Z.19 


6.28 


6.27 


4.97 


4.26 


ANTLIA 


79 


K on 

o.yu 


7/1 

o. ( 4 


o nfi 
Z.UO 


5.92 


5.90 


4.02 


2.58 


A2877 


65 


O.UU 


O 1 7 
O.l ( 


1 7C^ 

1.(0 


5.02 


4.99 


2.54 


2.66 


3C129 


70 


4. / 1 


o on 

z.yy 


1.00 


4.73 


4.70 


2.42 


1.02 


NGC5044 


35 


o.oy 


o oo 
z.ze 


1 Of^ 

l.zO 


3.61 


3.59 


2.35 


1.67 


NGC5846 


92 


o.o4 


O 1 o 

z.lz 


1 17 
1.1 < 


3.36 


3.35 


2.40 


3.40 


NGC1550 


15 


Z.OO 


1 f; 1 
l.bl 


n on 


2.54 


2.53 


1.36 


1.15 


A3571 


39 


z.z4 


1 y1 O 

1.4z 


n 70 
U. /8 


2.25 


2.22 


0.67 


0.45 


TRIANGUL 


95 


z.lo 


1 07 
l.O < 


n 7fi 
U. ( 


2.17 


2.13 


0.56 


0.42 


A2199 


51 


z.Uo 


1 on 
l.oU 


n 70 
U. ( z 


2.06 


2.03 


0.68 


0.50 


A2634 


60 


o no 
z.Uz 


1 oo 
l.z8 


n 71 
U. ( 1 


2.03 


2.01 


0.66 


1.03 


A3266 


17 


z.UU 


1 07 
l.z / 


n 7n 


2.01 


1.97 


0.46 


0.53 


A0496 


18 


i.yo 


1 oo 
i.zz 


n fi7 
U.O ( 


1.94 


1.91 


0.60 


0.49 


A0262 


5 


i.yz 


1 OI 

l.zl 


n fi7 
U.O / 


1.92 


1.91 


0.84 


0.64 


A0754 


23 


1.88 


1 in 

i.iy 


U.OO 


1.89 


1.85 


0.45 


0.66 


A0119 


2 


1.8 / 


1 in 
1.19 


U.OO 


1.88 


1.84 


0.49 


0.68 


IIIZw54 


12 


1 70 
1. /Z 


1 nn 


n cn 
U.OU 


1.73 


1.71 


0.53 


1.03 


A2319 


98 


l.DZ 


1 no 
l.Uo 


n c^7 
U.O / 


1.62 


1.59 
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11 
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0.24 
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NGC507 
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0.99 
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1.28 
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0.95 


0.94 
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Notes: (1) The HIFLUGCS cluster ID. (2) In units of 10""ph cm^^ s^V (3) In units of lO^^^^phcm^^ s^V 
(4) In units of 10~^^ergcm~^ s~^. (5) The flux within 0.1°, smoothed with a point spread function of 0.1°. 
(6) The flux within 1.0°, smoothed with a point spread function of 0.1°. (7) Signal to noise (S/N) within 7-200 
above 1 GeV. We estimate the signal to noise through V^obs ^off F-y/^/F-y + F^^n, where tobs = 3yrs is the 
observation time and A^s = 7000 cm^ is the effective area of Fermi-Lat at 1 GeV. 
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TABLE IX. Gamma-ray LP-IC-CMB and LP-BS flux within r2oo from brightest 50 clusters in HIFLUGCS catalogue. 
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82.88 



Notes: (1) The HIFLUGCS cluster ID. (2) In units of lO'^ph cm'^ s'^ (3) In units of 10~''phcm-2 8"^ (4) In 
units of 10~"'^^phcm~^ s~^. (5) In units of 10~^ergcm~'^ s~^. (6) The flux within 0.1°, smoothed with a point 
spread function of 0.1°. (7) The flux within 1.0°, smoothed with a point spread function of 0.1°. (8) Signal 
to noise (S/N) within r2oo above 1 GeV. We estimate the signal to noise through Viobs AcS F-^j ~J F~f + Ft^cr, 
where fobs ~ 3 yrs is the observation time and A^fi = 7000 cm^ is the effective area of Fermi-Lat at 1 GeV. 
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TABLE X. Galaxy cluster parameters for to all clusters in the HIFLUGCS catalogue. 
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,06 


0, 


,14 


A3827 


102 


451.10 


2 


.53 


0, 


.32 


18.70 


0, 


.31 


423 


0, 


.99 


0, 


,00 





0, 


,00 


0, 


,04 


0, 


,15 


A3888 


103 


720.64 


3 


.04 


0, 


.24 


32.40 


0, 


.61 


285 


0, 


.93 


0, 


,00 





0, 


,00 


0, 


,02 


0, 


,11 


A3921 


104 


429.52 


2 


.00 


0, 


.27 


9.20 


0, 


.40 


234 


0, 


.76 


0, 


,00 





0, 


,00 


0, 


,03 


0, 


,12 


HCG94 


105 


184.32 


1 


.49 


0, 


.46 


3.80 


0, 


.80 


42 


0, 


.53 


0, 


,19 


142 


0, 


,58 


0, 


,06 


0, 


,21 


RXJ2344 


106 


356.88 


2 


.17 


0, 


.35 


11.80 


0, 


.51 


91 


0, 


.72 


0, 


,32 


285 


0, 


,92 


0, 


.03 


0, 


,16 



Notes: (1) The HIFLUGCS cluster ID and luminosity distance (Dium) are taken from [H]. (2) The 
virial mass (M200) for each cluster is derived from the M500 mass in [35l [. We solve for M200 using 
M200 = M500 200/500 (c2oo(M2oo)/c50o(M2oo))^ [13. The virial radius (r2oo) is derived from M2oo- (3) 
The electron number density profile of each cluster follows either a single- or double-beta profile; 7ie(r) = 

y'X^i '^0,1(0)^(1 + where the central electron density, nc(0), core radius, rc, and slope, /3, are de- 

rived from [3^. (4) Half light radius (rhir.ca) of the flux from CR-7r°. (5) Half light radius (rhir,DM) of the 
flux from the BM-K' continuum and LP-IC-CMB emission. The flux has been boosted by substructures and 
includes the Sommerfeld enhancement for the LP model. 



